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Abstract 

We discuss Higgs boson decays in the CP-violating MSSM, and examine their phe- 
nomenological impact using cross section hmits from the LEP Higgs searches. This 
includes a discussion of the full 1-loop results for the partial decay widths of neu- 
tral Higgs bosons into lighter neutral Higgs bosons {ha — )• /ife/ic) and of neutral Higgs 
bosons into fermions {ha — )• //). In calculating the genuine vertex corrections, we take 
into account the full spectrum of supersymmetric particles and all complex phases of 
the supersymmetric parameters. These genuine vertex corrections are supplemented 
with Higgs propagator corrections incorporating the full one-loop and the dominant 
two-loop contributions, and we illustrate a method of consistently treating diagrams 
involving mixing with Goldstone and Z bosons. In particular, the genuine vertex cor- 
rections to the process ha — )• hbhc are found to be very large and, where this process is 
kinematically allowed, can have a significant effect on the regions of the CPX bench- 
mark scenario which can be excluded by the results of the Higgs searches at LEP. 
However, there remains an unexcluded region of CPX parameter space at a lightest 
neutral Higgs boson mass of ~ 45 GeV. In the analysis, we pay particular attention 
to the conversion between parameters defined in different renormalisation schemes and 
are therefore able to make a comparison to the results found using renormalisation 
group improved/effective potential calculations. 
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1 Introduction 



High energy colliders, at present and in the past, have given the search for Higgs bosons a 
high priority. The LEP and Tevatron experiments, in particular, have been able to turn the 
non-observation of Higgs bosons into constraints on the Higgs sector, which have been very 
useful in reducing the available parameter space of some of the most popular particle physics 
models, such as the Standard Model (SM) [1] and the Minimal Supersymmetric Standard 
Model (MSSM) [2]. For first results on the Higgs searches at the LHC, see Refs. [3,4]. 

However, MSSM scenarios involving CP violation in the Higgs sector, which induces a 
mixing of all three neutral Higgs bosons, can prove particularly difficult to restrict using the 
Higgs search data. This is due to the fact that the CP violation can result in suppressed 
couplings of the lightest Higgs boson to two gauge bosons and to the non-standard decay 
mode of a heavier SM-like Higgs boson into a pair of light Higgs bosons, resulting in an 
experimentally rather challenging final state. The CPX benchmark scenario [5] is an example 
of such a situation in the MSSM. In the original combined LEP analysis by the LEP Higgs 
Working group and the LEP collaborations (LHWG), it was found that substantial regions 
of the CPX parameter space could not be excluded [2] where the lightest Higgs mass is 
substantially below the limit on the Standard Model Higgs mass [1] of Mh = 114.4 GeV. 

In this paper, we will present complete one-loop results for the decay widths of neutral 
Higgs bosons into lighter neutral Higgs bosons (Higgs cascade decays) and the decay widths 
of neutral Higgs bosons into fermions in the CP- violating MSSM. The results are obtained in 
the Feynman-diagrammatic approach, taking into account the full dependence on the spec- 
trum of supersymmetric particles and all complex phases of the supersymmetric parameters. 
The genuine vertex contributions are supplemented with two-loop propagator-type correc- 
tions, yielding the currently most precise prediction for this class of processes. One-loop 
propagator-type mixing between neutral Higgs bosons and Goldstone and Z bosons is also 
consistently taken into account. 

Both of these calculations require loop corrections to the neutral Higgs mass matrix M, 
which are well known for the real and complex MSSM and are frequently used to add propaga- 
tor corrections to processes involving external neutral Higgs particles. These corrections are 
incorporated in the two main public codes for calculating the complex MSSM Higgs sector, 
FeynHiggs [6-10] and CPsuperH [11,12]. FeynHiggs is based on the Feynman-diagrammatic 
approach and on-shell mass renormalisation while CPsuperH is based on a renormalisation 
group improved effective potential calculation and DR renormalisation. Therefore, to com- 
pare between these results it is necessary to perform a parameter conversion. We shall discuss 
this issue in Sect. 5. We also investigate the numerical impact of parametrising the neutral 
Higgs self-energies (in the Feynman-diagrammatic approach) in terms of the MS top mass, 
rather than the on-shell top mass, which is formally a 3-loop effect. 

The Higgs cascade decays often dominate the Higgs decay width where they are kinemat- 
ically allowed. They directly involve the Higgs self-couplings, the observation and measure- 
ment of which is a crucial goal for the experimental confirmation of the Higgs mechanism. We 
will present two momentum-dependent approximations for the loop-corrected triple Higgs 
couplings, which can be used, for instance, for predictions of the Higgs production process 
e+e- ^ Zhaha at the ILC [13] or CLIC [14]. 

The genuine vertex corrections to the triple Higgs decay can be very large. In the 



1 



MSSM with real parameters, the leading Yukawa vertex corrections and the complete 1- 
loop vertex corrections have been calculated [15-21]. However, for the complex MSSM, 
previous to our result, first described in Ref. [22], only effective coupling approximations 
were available [23,24], as provided by the program CPsuperH [11,12]. The genuine vertex 
corrections we present will be incorporated into the code FeynHiggs. As we will demonstrate, 
the /i2 — hihi decay width has a critical influence on the size and shape of one of the regions 
of CPX parameter space which the LEP Higgs search results are unable to exclude. 

The fermionic decay modes of the neutral Higgs bosons are crucially important to collider 
phenomenology. These modes have been used when obtaining a lower bound on the Standard 
Model Higgs mass [1] and to exclude significant regions of the MSSM parameter space 
[2,25,26]. In particular, an accurate prediction for the Higgs decay to b-quarks has been 
vital for these analyses, since, for Standard Model Higgs bosons with mass less than about 
130 GeV and for most SUSY scenarios, ha — ?■ hh is the dominant decay mode. The decay 
to r-leptons can also be very important for Higgs searches, as demonstrated for instance for 
various benchmark MSSM scenarios in the high tan/3 region at the Tevatron [27]. 

In the Standard Model, the fermionic decay width is extremely well known (for a review, 
see e.g. Ref. [28] and references therein), and the treatment of higher-order QCD (gluon- 
exchange) and QED corrections can be taken over to the MSSM case. The SUSY QCD 
corrections can be sizable for the ha — ?■ hh decay and should be resummed (see, for example, 
Ref. [29], for an investigation into these effects). Results supplemented with leading 2-loop 
propagator corrections [30] and full electroweak contributions [31] are also available in the 
MSSM with real parameters. 

Predictions for the ha — > // decay widths for the Standard Model and the MSSM with 
real parameters can be obtained from the programs HDECAY [32] and HFOLD [33]. For 
the complex MSSM, the program CPsuperH [11,12] is available. It is based on calculations 
involving effective haff couplings, as described in Ref. [24]. 

The program FeynHiggs [6-10] calculates the ha — > // decay width using the Feynman- 
diagrammatic approach, including the most significant QCD corrections, resummed SUSY 
QCD corrections and propagator corrections incorporating the full neutral Higgs self-energies. 
This calculation is valid in the real and complex MSSM. The full 1-loop electroweak vertex 
corrections presented here have recently been incorporated into FeynHiggs. 

The corrections to the Higgsstrahlung and Higgs pair production processes at LEP in the 
MSSM with real parameters have been studied in Refs. [34-39] and the CP-violating MSSM 
in Refs. [40-44]. In the present paper, we will investigate the t,t,h,h corrections to these 
production processes in the Feynman-diagrammatic approach in the CP-violating MSSM, 
and supplement these with full propagator-type corrections. This type of corrections were 
not included in the Feynman-diagrammatic analysis of the CPX scenario in Ref. [2]. 

The parameter region in the MSSM with complex parameters that could not be excluded 
with the Higgs searches at LEP, characterised by a rather light Higgs boson with a mass 
of about 45 GeV and moderate values of tan/3, persists also in view of the present search 
limits from the Tevatron [27]. This parameter region will be difficult to cover also with the 
standard Higgs search channels at the LHC [45-47], while it can be thoroughly investigated 
at the ILC [13]. The phenomenology of scenarios with such a light Higgs boson has recently 
found considerable interest in the literature, see Refs. [48-56] for discussions of other (non- 
standard) possible LHC search channels to access this parameter region. 
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In the present paper we make use of our improved theoretical predictions for the Higgs 
branching ratios into a pair of hghter Higgs bosons and into a fermion pair to examine their 
impact on the parameter region with a hght Higgs boson left unexcluded by the LEP Higgs 
searches. For this purpose we employ the topological cross section limits obtained at LEP, 
as implemented in the program HiggsBounds [57,58]. We investigate the sensitivity of the 
excluded parameter region with respect to variations in the parameters of the CPX scenario. 
This analysis updates and considerably extends our previous results reported in Ref. [22]. 
We then compare our results to the results obtained with the code CPsuperH, using various 
ways of performing the parameter conversion. 

The paper is organized as follows: After introducing complex parameters in Sect. 2 and 
the CPX scenario in Sect. 3 we discuss contributions to the Higgs masses and mixings includ- 
ing also resummed SUSY QCD corrections in Sect. 4. In Sect. 5 we focus on the conversion 
between different renormalization schemes as well as on the effect of a different parameteri- 
zation of the top quark mass. In Sect. 6 and in Sect. 7 we discuss the Higgs cascade decay 
and the Higgs decay into SM fermions, respectively, and the different contributions to their 
partial decay widths and possible approximations. After the investigation of the partial de- 
cay widths we turn our focus particularly on the branching ratios of the Higgs cascade decay 
processes in Sect. 8. In Sect. 9 Higgs production channels which were relevant at LEP are 
investigated. Finally, in Sect. 10 we discuss the phenomenological impact of the improved 
predictions obtained in this paper. We investigate in particular the parameter dependence 
of the CPX scenario and we perform a thorough comparison with the results obtained with 
the program CPsuperH. Sect. 11 contains our conclusions. 

2 The MSSM with complex parameters at tree level 

In its general form, the MSSM allows various parameters to be complex. This includes the 
trilinear couplings Af, the Higgsino mass parameter /i, the gluino mass parameter M3 and 
the soft SUSY breaking parameters Mi and M2 from the neutralino/chargino sector. These 
complex parameters can induce CP violation. Below we list the relevant quantities to fix 
our notation, which closely follows that in Ref. [7]. 
We write the two MSSM Higgs doublets as 



where vi and V2 are the vacuum expectation values, and tan/3 = V2/V1. Here we have made 
use of the fact that the MSSM Higgs sector is CP-conserving at lowest order, i.e. complex 
phases occurring in the Higgs potential can be rotated away (or vanish via the minimisation 
of the Higgs potential). 

The tree level neutral mass eigenstates h, H, A, G are related to the tree level neutral 
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fields 01, Xi) X2 through a unitary matrix 
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cos a 
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sin a 








02 








— sin /3n 


COS /3n 


Xi 








cos /3n 


sin /3n / 


\X2j 



(2) 



in which the CP-even eigenstates 0i,02 do not mix with the CP-odd eigenstates Xi)X2- 
Unless otherwise stated, h, H, A, G will always represent tree level neutral (mass eigenstate) 
fields throughout this paper. At tree level, the off-diagonal mass terms must vanish, leading 
to the condition (3^ = (3. 

The Higgs sector at lowest order is given in terms of two independent parameters (besides 
the gauge couplings), conventionally chosen as tan/3 and either rriA or mH±. Since CP 
violation can be induced via potentially large higher-order corrections, in general all three 
neutral Higgs bosons will mix once higher-order corrections are included, so that the CP-odd 
A boson is no longer a mass eigenstate. For the general case of the MSSM with complex 
parameters it is therefore convenient to use mH± as input parameter. In our notation lower- 
case Higgs masses indicate tree-level masses, while upper case masses refer to loop-corrected 
masses. 

We write the squark mass matrices as 



Ml 



m: 



+ M| cos 2/3(/|-Qgs2 



iTigX* 



nig Xq 



M^^ + m' + M'z cos 2 (3QgSt 



where 



Aq — /i*{cot (3, tan/3}. 



(3) 



(4) 



and cot /3 or tan/3 applies to u-type or d-type quarks, respectively. The eigenvalues of eq. (3) 
are 

1 



91,2 



2 L 



Ml + Ml 



/|M|cos2/3 



T J [Ml - Ml + M| cos2/3(/| - 2g,4)]2 + 4m2|X, 



(5) 



In the complex MSSM, the trilinear coupling Aq and the higgsino mass parameter /i can 
have non-zero complex phases. The mass matrix Mg can be diagonalised by the matrix Ug. 
Here 



Q2 



where 



(6) 



and Cq is real, Sq is complex, and c? + = 1. 

The coefficient of the gluino mass term in the Lagrangian, M3, is in general complex. 
The gluino mass is given by nig = \Ms\, while the phase 0m3 can be absorbed into the gluino 
fields [59]. The phase of M3 thus appears in the quark-squark-gluino couplings. 

For the chargino mass matrix we use 



chargino 



M2 V2 sin (3 Mw 

cos /3 Mw fJ- 



(7) 
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which includes the soft SUSY-breaking term M2, which can be complex. For the neutralino 
mass matrix we use 




— M^SwCos/3 Mz sin l3\ 
Mz Cw cos l3 Mz Cw sin /3 



■neutralino 



-M2 



-/i 
-/i / 



8 



which includes furthermore the soft SUSY-breaking term Mi, which can also be complex. 

It should be noted that not all phases mentioned above are physical, but only certain 
combinations. In particular, the phase of the parameter M2 (chosen by convention) and, as 
mentioned above, the phase appearing in the Higgs sector at lowest order can be rotated 
away. 

3 Phenomenology and the CPX scenario 

CP-violating effects, which can enter the Higgs sector via potentially large higher-order 
corrections, can give rise to important phenomenological consequences. CP phases in the loop 
corrections to the Higgs particles can have a large impact on the predictions for the masses 
(all three neutral Higgs bosons mix in the CP- violating case) and the Higgs couplings [7,41, 



Studies of the possible impact of CP-violating effects on the MSSM Higgs sector have 
often been carried out in the CPX benchmark scenario [5]. As input values for the CPX 
scenario we use in this paper 

• mt = 173.1 GeV 

• MsusY = 500 GeV(= M^""*^^"^^ = M?;;-^^'"") 

• /i = 2000 GeV 

• iM-sl = 1000 GeV 

• M2 = 200 GeV, Ml = f^Ms (see e.g. Ref. [63]). 

• lA^'^-^hciii = i^^i = GeV 

• </)^on-shell = (pA,, = 4>M-i = f 

• Mh± < 1000 GeV 

With the phases of the parameters At,b and M3 set to the maximal value of 7r/2 and the 
relatively large value of /i, this scenario has been devised to illustrate the possible importance 
of CP-violating effects. 

The above values differ from the ones defined in Ref. [5] in the following ways: Firstly, 
we use an on-shell value for the absolute value of the trilinear coupling At and the soft 
SUSY breaking mass parameters M^ and M^^, rather than DR values, and we therefore 
use a numerical value of \At\ that is somewhat shifted compared to that specified in Ref. [5] 



60-62]. 
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in order to remain in an area of parameter space with similar phenomenology (the value 
specified in Ref. [5] is \Af^\ = 1000 GeV). Secondly, we use rrit = 173.1 GeV, which was 
the world average top-quark mass in March 2009 [64]. 

We use an on-shell definition of At, and Mg^ since this is the natural choice for a 
Feynman-diagrammatic calculation. We will discuss how to convert between the different 
parameter definitions in Sect. 5. For the purposes of this discussion, we use a second scenario 
using the parameter values given above, except with At, Ml, Mg^ defined according to the 
DR scheme at the scale Ms := ^/M|^^y + and with \Af^{Ms)\ = 1000 GeV(= \Ab\), 
(l)^m{Ms) = f , M^{Ms) = 500 GeV and MPf (M^) = 500 GeV, which we will call the 

CPX scenario (i.e. this scenario is more similar to that in Ref. [5]). 

The LEP Higgs Working Group study [2] of the CPX scenario also used \Af^\ = 
1000 GeV. The majority of its analyses were performed using rrit = 174.3 GeV. We will 
investigate the dependence of our results on rrit in Sect. 10.2. 

It should be noted that there are constraints on the CP phases in the complex MSSM 
from experimentally measured upper limits on electric dipole moments, such as those of the 
electron and neutron (for a recent discussion, see e.g. Ref. [65]). These provide particularly 
significant constraints on the CP phases in the first two generations. The constraints on the 
phases of the third generation are less restrictive. In the definition of the CPX benchmark 
scenario, existing bounds on CP phases were taken into account, see Ref. [41] for more 
details. 

4 Loop corrections to the Higgs masses and Higgs mix- 
ing matrices 

Higher-order corrections to Higgs masses and mixing properties are known to be very im- 
portant for the phenomenology of the MSSM Higgs sector, see Refs. [66-68] for reviews. 

In the MSSM with real parameters, the full 1-loop result [31, 34, 69-74] and the dom- 
inant 2-loop corrections [9, 75-89] have been calculated, and the tan /3-enhanced terms 
0(Q;fe(as tan/3)") have been resummed [29,90-94]. A full 2-loop effective potential cal- 
culation is known [95-102]^. In addition, some dominant 3-loop contributions have been 
calculated [103-105]. 

In the complex MSSM, 1-loop corrections from the fermion/sfermion sector and some 
leading logarithmic corrections from the gaugino sector and the dominant 2-loop results have 
been calculated in the renormalisation group improved effective potential approach [41,60,62, 
97,100,104,106-108]. In the Feynman-diagrammatic approach, leading 1-loop contributions 
have been obtained in Ref. [61, 109], and the full 1-loop result has been calculated in Ref. [7]. 
At 2-loop order, the 0{atas) corrections are available [110]. 

Most of these results for the complex MSSM have been incorporated either into the 
public code FeynHiggs [7-10,111-113], which uses the Feynman-diagrammatic approach, or 
the public code CPsuperH [11,12], which uses the renormalisation group improved effective 
potential approach^. 

^In principle, the effective potential calculation is also applicable to the complex MSSM. 

^Unless explicitly stated otherwise, 'FeynHiggs' will refer to FeynHiggs version 2.6.5 and 'CPsuperH to 
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In this paper, when calculating the Higgs masses and mixings, we will use renormalised 
neutral Higgs self-energies calculated by FeynHiggs, to take advantage of the fact that it 
includes the complete 1-loop result and 0{atas) corrections of Ref. [110] with full phase 
dependence. FeynHiggs additionally allows the option of including sub-leading 2-loop cor- 
rections which are known so far only for the MSSM with real parameters [76,82,83,85,86]. 
If the user wishes to apply these corrections in an MSSM calculation with complex phases, 
FeynHiggs evaluates these corrections at a phase of and vr for each complex parameter, 
then an interpolation is performed to arrive at an approximation to these corrections for 
arbitrary complex phases. However, this prescription can be problematic in a rather 'ex- 
treme' scenario like the CPX scenario. In fact, it can happen in this case that one of the 
combinations of real parameters needed as input for the interpolation turns out to be in an 
unstable region of the parameter space where the reliability of the perturbative predictions 
is questionable. This would skew the interpolation towards the unstable values. Therefore, 
unless otherwise stated, we will use in the present paper the leading 0{atas) corrections 
to the Higgs self-energies from FeynHiggs, but not the sub-leading 2-loop corrections. A 
discussion of the incorporation of the subleading 2-loop contributions via the interpolation 
from the results for real parameters is given in Sect. 10.3. As discussed in more detail below, 
besides the irreducible 2-loop contributions of 0{atas) we do incorporate into our results 
higher order tan /3-enhanced terms (for arbitrary complex parameters), which we take into 
account by introducing an effective b-quark mass. 



4.1 Determination of neutral Higgs masses 

In general, the neutral Higgs masses are obtained from the real parts of the complex poles of 
the propagator matrix. In the determination of the Higgs masses, we neglect mixing with the 
Goldstone and Z bosons as these are sub-leading 2-loop contributions to the Higgs masses. 
We therefore use a 3 x 3 propagator matrix A(p^) in the {h, H, A) basis. 

In order to determine the neutral Higgs masses we must first find the three solutions to 

\A-\p')\=0, (9) 

which, in the case with non-zero mixing between all three neutral Higgs bosons, is equivalent 
to solving 

where i = h, H or A. The propagator matrix is related to the 3x3 matrix of the irreducible 
2-point vertex-functions F2{p'^) through the equation 



where 



[-A{p')y' = hip') = z [/II -M(/)] , (11) 



M{p') = I -f:,H{p') ml - tuHip') -±ha{p') I • (12) 
-thAip'') -tHAip") 




CP super H version 2.2 throughout this paper. 
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As before, mh,mH,mA refer to the tree level masses. Sjj(p^) are renormalised Higgs self- 
energies. The explicit form of the counterterms used in this paper is given in App. B. For 
the majority of these renormalisation conditions, we use the on-shell scheme. However, we 
shall use DR renormalisation for the Higgs fields (see Ref. [7]). If there is CP conservation, 
^/ia(p^) = ^HAip"^) = 0, and the CP-even Higgs bosons h,H do not mix with the CP-odd 
Higgs boson A. 

In general, the renormalised Higgs self-energies can be complex, due to absorptive parts. 
Therefore, the three poles of the propagator matrix Ai^ can be written as 

Ml = Ml-tM,^W^^, (13) 

where M^^ is real and is interpreted as the loop-corrected (i.e. physical) mass, Wh^ is the 
Higgs width, and a = 1, 2, 3. 

In the MSSM with complex parameters, the loop-corrected masses are labelled in size 
order such that 

Mh, < Mh, < M;,3. (14) 

In the CP-conserving case, the masses are labelled such that the CP-even loop-corrected 
Higgs bosons have masses and Mh with < and the CP-odd loop-corrected 
Higgs boson has mass Ma (the numerical value of Ma is affected by loop corrections only if 
Mh± is chosen as an independent input parameter). 

Solving eq. (9) with full momentum dependence involves an iterative procedure, since the 
self-energies also depend on the momentum. In order to deal with the complex momentum 
argument it is convenient to use an expansion about the real part of the pole, M|^, such 
that 

±,kiMl) ^ ±,,iMl)+tlm[Ml]±'^,iMl), (15) 

with j = h,H,A and k = h,H,A. We obtain Ejfc(M^J and f:'ji^{MlJ (where the prime 
indicates the derivative w.r.t. the external momentum squared) from FeynHiggs [6-10]. For 
each ha, we use a momentum-independent approximation to obtain an initial value for the 
iteration. This solution is then refined using 

Al2>+il = ath eigenvalue of M(7Wj^'''l), (16) 

where the eigenvalues have been sorted into ascending value, according to their real parts. 
We check the validity of the truncation of the expansion in eq. (15) by performing an iteration 
using an expansion up to second order: 

±,kiMl) ^ ±,kiMl) + tlm[Ml]^kiMl) 

+ l{i^^[^U^%iMl)y, (17) 

and confirming that the resultant Higgs masses show no significant changes. 
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4.2 Wave function normalisation factors 



In order to ensure that the S-matrix is correctly normahsed, the residues of the propagators 
have to be set to one. We achieve this by including finite wave function normalisation factors 
which are composed of the renormalised self-energies. These 'Z-factors' can be collected in 
to a matrix Z where 



lim 



lim 
lim 



i 



Z To Z 



Z ■ To ■ 



hh 



HH 



Ml 



z ■ r, ■ z 



AA 



(18) 
(19) 
(20) 



such that 



(21) 



where f is a one-particle irreducible n-point vertex-function which involves a single external 
Higgs ha, and ha, hf,, he = some combination of hi, /12, h^. 
The matrix Z is non-unitary. We write it as 



ZhZhh 

'ZAZAh 



ZaZah 




(22) 



We find the elements of Z by solving eq. (18), which gives 
1 



Zh 



a 



ZhH 
ZhA 



AhA 



A 



hh 



^Hh 



^HA 



Za = 

ZAh 

Zah 



A 



hA 



Aaa 
Aha 



i(p2) 



lie 



p'^=Ml 



A 



AA 



(24) 
(25) 



'he 



We choose ha = hi, hh = /i2 and he = h^. Zi is, in general, complex. Other choices for 
the Z-factors are possible, such as that in Ref. [114], where we use the limit = M\^^,^. 
However, this does not allow the same freedom for choosing a, b, c. 

Since the elements of Z involve evaluating self-energies at complex momenta, we use 
again the expansion given in eq. (15). In order to make sure that the neglected higher order 
terms in eq. (15) are small, we also calculate Z using eq. (17), and check that this does not 
significantly change the result. 

The wave function normalisation factors are included in the calculation by multiplying 
the irreducible vertex factor F by Z once for each external Higgs boson involved in the 
process. 
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Figure 1: The Goldstone and Z boson propagator corrections to the hi — // decay, where 
hi = h, H OT A 



4.3 Goldstone or gauge bosons mixing contributions to the Higgs 
propagators 

A complete 1-loop prediction for a process involving a neutral Higgs propagator in the MSSM 
with complex parameters will, in general, contain terms involving the self-energies S/iG, ^hg, 
Tjag and T^hz, ^hz, ^az, such as those shown in Fig. 1. These terms are required to ensure 
that the 1-loop result is gauge-parameter independent and free of unphysical poles. As we 
will illustrate for an example, it is essential to treat these mixing contributions strictly at 
one-loop level in order to ensure the cancellation of the unphysical contributions (some care 
is necessary to achieve this, since the loop corrected masses used for the external particles 
and the Z-factor prescription outlined above automatically incorporate leading higher-order 
contributions). 

As an example, we consider diagrams involving mixing contributions for a neutral Higgs 
decaying to two fermions, as in Fig. 1. We use here the lowest order Z-boson propagator 
with explicit gauge parameter dependence in the gauge, 

P/^pA i Pt^Pu i^z 



and the G-boson propagator 



(27) 



The vertex F^g^^™ involving on-shell fermions is related to FqJ^ by 

P.r^r = -^MzFg-. (28) 

The relation between the hG and hZ self-energies is given in eq. (130) below. Using this, 
we can express the h ff decay (Fig. 1 with hi = h) as 

■A "ptree 

'^^"-^f'^f^^'SS + 'P-^-^(p \.(;l""4fe) rSr = -rf (29) 

where denotes the momentum of the propagator involving the neutral Higgs boson. Note 
that the expression in eq. (29) does not contain a pole at = M'^^z- Similarly, eq. (132) 
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Figure 2: SUSY QCD corrections to the relation between the bottom quark mass and the 
bottom Yukawa couphng induced by gluino and sbottom quark loops, which can be enhanced 
at large tan/3, {x,y = 1,2) 



below gives the relation between the AG and AZ self-energies. When substituted into the 
expression for the decay A — )■ // via a self-energy (Fig. 1 with hi = A), this gives 



^^^^^^ V^-Mk/ ^"+^^ ^^^^^ \^ip^-MUzf -^-^ 



■ptree / i\/f2 f 

Gff / V /^2\ / 2 2 \ r / 2\ ^^^Z^Z 



p2 y^Ao{p')-{p'~Ofo{p')^^^J. (30) 

and the quantity /o is defined in eq. (128) below. The expression above shows that it is 
essential to use the tree-level mass for the incoming momentum, i.e. p"^ = m\, in order to 
ensure the cancellation of the unphysical pole at p"^ = M'^^z- Therefore, in the following, 
we treat the contributions involving mixing between hi and G, Z bosons strictly at one-loop 
order, which implies, in particular, evaluating those contributions at an incoming momentum 
corresponding to the tree level mass, rather than the loop corrected mass. 

4.4 Resummation of SUSY QCD contributions 
4.4.1 The Am;, correction 

The tree level relation (m^ = XbVi) between the bottom quark mass and the bottom Yukawa 
coupling Afe receives large tan /3-enhanced radiative corrections, which need to be properly 
taken into account [29,90-94,115-117]. In SUSY QCD, such contributions arise from loops 
containing gluinos and sbottoms, as shown in Fig. 2. For heavy SUSY mass scales the 
interaction of the neutral Higgs bosons with bottom quarks can be expressed in terms of an 
effective Lagrangian [29] 



-nil H 7 -^22 

*/3 



h.c, (31) 



where the shorthand t/3 = tan/3 has been used. Accordingly, the relation between the bottom 
quark mass and the bottom Yukawa coupling receives the loop-induced contribution Am^ 



11 



such that 



"rrib = XbVi (1 + Amfc) 



(32) 



We consider here the general case, in which Am^ is allowed to be complex. Inserting the 
relation (eq. (32)) and neglecting the terms involving Goldstone boson contributions leads 
to 



cfT 



+ 



+ 



t 



y + ^753; 



^hbb^ 



1 + —y + i-f^x 1 - — 

*/3 V */3/J 



1 



where ta = tana. The quantities x,y are real and given by 

ImAmf, 
1 + ReArrtft ' 
ReAmft + xlmAmfe, 




^l5b^ 1 b+..., 



(33) 



X 



y 



and v'^^^, v'^g, vfg 



tree 



are defined by 



A 



(0) 



^/2 



with 



A 



(0) 



nib/vi 



A 



(0) 



V2 



V2 



{il^){-sp)A 



b + ... 



rribe 



(34) 

(35) 
(36) 

(37) 



Note that, in this convention, u^^^ contains a 75 dependence. 

In order to find Am;,, we perform a Feynman-diagrammatic calculation of the leading 
1-loop gluino contributions to the hi — )■ bb decays, using the p"^ = approximation and 
i = h,H, A. Comparing this calculation for the renormalised decay width to the 1-loop 
expansion of eq. (33) yields for the contribution of gluino and sbottom loops to Amb 



Ami 



I (a, b, c) 



,7717 , m 



abLog (^) + acLog (f ) + 6cLog (f ) 
(a — c) (c — b) {b — a) 



(38) 
(39) 



In the hi — )■ bb decay, diagrams involving charged higgsinos also contain tan /3 enhanced 
contributions [29,118]. We treat these analogously to the Am^ corrections above. Compar- 
ison with the 1-loop Feynman-diagrammatic calculation in the complex MSSM leads to 
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Figure 3: Predictions for the neutral Higgs masses, M/i^, M^^, Mh^, obtained using the 
effective bottom quark mass employed for the numerical analyses in this paper, (iii) off = 

m^^''^^(m^'^)/|l + Am™ I in comparison with two other choices for the bottom mass, (i) 
^MS,SAfj^^^-j and (ii) n\^'^^\m'^^). The results are shown for the CPX scenario with 
Mh± = 140 GeV. 



where 



x(0),2 

A?^ = !!^ = ^^ . (42) 



V2 V2sf}SwMw 

The effective Lagrangian of eq. (31) properly resums the leading tan /3-enhanced gluino 
and higgsino contributions given above [29]. In our calculation we therefore use a Am^ 
correction of 

Amb = Ami + ^"^fe (43) 
= l^f^*M;t,I {ml,ml,ml) + ^A^^^HpI (m|,m|, \^^\') . (44) 

It is also possible to incorporate effects from loops involving winos into Arrif, as in Ref. [29] 
(or even winos and binos as in Ref. [118]). We do not include these, since they are numerically 
small [29] and, in the CPX scenario, are less important than the higgsino contributions. Since 
we will explicitly calculate the 1-loop diagrams involving winos and binos when calculating 
the full 1-loop ha — bb decay width (the Arub contributions are subtracted at one-loop order 
such that a double-counting from the Arub resummation is avoided), the effect of leaving 
them out of the Arrib contribution is of sub-leading 2-loop order. For the scale of as in 
eq. (44) we choose the top-quark mass, i.e. we use as^mf) in Anib- 

4.4.2 The use of an effective 6-quark mass in the calculation of the neutral 
Higgs self-energies 

It is desirable to incorporate leading tan /3-enhanced contributions also into the calculation 
of the Higgs self-energies. A direct application of the effective Lagrangian of eq. (31) within 
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loop calculations is in general not possible, as it would spoil (among other things) the UV- 
finiteness of the theory. It is therefore convenient to absorb the leading tan /3-enhanced 
contributions into an effective bottom quark mass that is used everywhere in the calcula- 
tion, although this procedure does not reproduce the decoupling properties of the effective 
Lagrangian of eq. (31) in the limit where Mh± 3> Mz- It has been shown in Refs. [76,85] 
that the one-loop result with an appropriately chosen effective bottom quark mass in general 
approximates very well the result containing the diagrammatic two-loop contributions. 

Because of the large value of /i in the CPX scenario, contributions from the bottom 
/ sbottom sector can be important already for moderate values of tan/3. We use in the 
following an effective bottom quark mass that is defined as 

DR,5M/ oS\ 

mb,cs = FHI ' ^5 

|1 + Am^^l 



where 'DR' or 'OS" indicates the renormalisation scheme in which the mass is defined i.e. 
the DR or 05" renormalisation scheme respectively, and 'SM' indicates that only the Stan- 
dard Model contribution is included. Am™ is the Ami, correction calculated internally by 
FeynHiggs for use in its Higgs decays and Higgs production cross sections. 

Since in the CPX scenario |1 -|- Amb\ > 1, the incorporation of the tan /3-enhanced Am^ 
corrections leads to a reduction of the numerical value of mt^eff in this scenario and thus 
to numerically more stable results. This is illustrated in Fig. 3, where the results for the 
neutral Higgs masses obtained using the effective bottom quark mass defined in eq. (45) 
are compared with the predictions arising from choosing m^^''^*^(mb) or m^^'^^"^ (mf^) for 
the bottom quark mass. For large values of tan/3, depending on the choice of the bottom 
mass, an onset of very large corrections from the bottom / sbottom sector is visible in 
the predictions for the Higgs masses. Since perturbative predictions are not reliable in the 
parameter region where these large corrections from the bottom / sbottom sector occur, we 
limit the numerical analyses in this paper to the region tan /3 < 30. (Without resummation 
of the tan /3-enhanced contributions as in eq. (45), the onset of very large corrections could 
already occur at lower values of tan /3.) 



5 Conversion of parameters between on-shell and DR 
renormalisation schemes 

As mentioned above, higher-order contributions in the Higgs sector of the MSSM have been 
obtained using different approaches. The results implemented in the public code Feyn- 
Higgs [6-10] are based on the Feynman-diagrammatic approach employing the on-shell renor- 
malisation scheme, while the results implemented in the public code CPsuperH [11,12] are 
based on a renormalisation group improved effective potential calculation employing DR 
renormalisation. As the parameters in the two renormalisation schemes are defined differ- 
ently, a parameter conversion is necessary for a meaningful comparison of results obtained 
in the two schemes. 
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Figure 4: The diagrams used to calculate the shifts shown in eqs. (48) - (53), (which convert 
between DR and on-shell parameters) at 0{as). (x = 1, 2, ?/ = 1, 2, z = 1,2) 

5.1 Parameter shifts 

Since both schemes incorporate partial 2-loop contributions, a parameter conversion of the 
top/stop sector parameters (that enter at the 1-loop level), is required. In Refs. [77,82] 
this issue has been discussed for the case where all the MSSM parameters are real. In the 
following we consider the general case of arbitrary complex parameters. We can obtain the 
leading terms at 0{as) from considering loops involving gluons, gluinos, stops and tops as 
shown in Fig. 4. In order to obtain the leading terms at 0{at), we must also consider loops 
involving neutralinos, charginos, Higgs bosons, Goldstone bosons, sbottoms and b-quarks, 
as shown in Fig. 5. 



We label the difference between the parameters p in the different renormalisation schemes 
by Ap, where 

pDR ^ p°--^-'' + Ap{fi,,^). (46) 

Since the DR parameters depend on the renormalisation scale /iren, the shift Ap is also a 
function of /^ren- The parameter shift Ap is related to the counterterms by 

Apiflrcn) = - 6p^ = ^p^'^-^hcU _ ^^pon-shclljdiv ^ ^^^^ 

where the superscript 'div' denotes that only terms proportional to — 7_e + log(47r) 

are kept. Therefore this means that Ap(/ii.cn) = [5^°""'^'^'^"]'^'^, where the superscript 'fin' 
denotes the finite pieces remaining once terms proportional to — 7s + log(47r) have been 
subtracted out. 

For the stop sector, we can directly adapt the counterterms used in FeynHiggs in Ref. [110] 
(see App. B for the explicit form of the counterterms) to get the parameter shifts 

AMI = -2mtAmt + f/*if/nAm| + U;jJ2iAm\ + Ul-^U2iAYi + U^^UiiAY-* , (48) 
AM^ = -2mtAmt + Ul^U^2Aml + t/*2?722Am| + Ul^U22AYi + U;^U^2AY~* , (49) 



15 



hi,G H ,G Xn 

t t t t 



t t \^ t \.^y t V.,,/ 

t b f h 

l^x Ox 

hi,G ,tz,bz 

t ^ b . 

-^"O^" --~-cy- 

tx ty tx ^ tx ty tx 2*^^ 

XU -y — 

n ■^z 



hi,G _ H ,G _ 



Figure 5: The diagrams used to calculate the shifts shown in eqs. (48) - (53) (which convert 

between DR and on-shell parameters) at 0{at). {x = 1,2, y = 1,2, z = 1,2, i = 1, 2, 3, n = 
1,2,3,4) 

AA* = e-^'^^*(A|A|-2|A|A0Aj, (50) 
AKt = - (A* - /i cot (3) Anit + nit cot /3A/i 

+[/i*if/i2 Am| + [/*it/22Am| + f/^if/asA^ + f/2*if/i2 AF/, (51) 

A|A| = —ReW'^^'^AKt], (52) 

= \-.\m[e'^^^AKt], (53) 



where 



Amt = ^Re[m(Sf(m2) + Sf(m2)) + Sj(m2) + S[(m2)]''", (54) 

Ami = (55) 

Am| = Re[E,Jml)f\ (56) 

1 ^ 1 fin 



AF. 



2 ^ReS,-^,(m|) + ReS,^,(m|)J , (57) 

and the components of the top self-energy are defined by T,{p) = ;^uj^T,^{p'^) +|/u;+S^(p^) + 
a;_E'(p^) + u^Ti^Ip^), and Re indicates that the imaginary parts of the loop integrals are 
discarded. 

In the following we will define the DR parameters at the scale ^^en = ^/ ^Iusy + ~- 
Ms, see Ref. [11] and the discussion in Ref. [77]. We evaluate the strong coupling constant 
at the scale of the top mass, ^^(m^). Since the two-loop corrections of 0{atas) and 0{af) 
implemented in FeynHiggs have been obtained using a Yukawa approximation (the corre- 
sponding contributions are implemented in CP super H oiAy to leading logarithmic accuracy). 



16 



we employ the same kind of approximation to derive the parameter conversions at 0{as) 
and 0{at). In particular, we neglect the D-terms in the stop mass matrix (i.e., the terms 
proportional to M|), the D-terms and the b-quark mass in the sbottom mass matrix, and 
we neglect terms proportional to Mi, M2, Mz, Mi^ in the neutralino and chargino mass 
matrices. In addition, we make the approximation Mh± = when deriving the parameter 
shifts. 

For the evaluation of the parameter shifts at 0{at) we also need to consider a shift to the 
Higgsino mass parameter /i, since FeynHiggs takes fi^^{mt) as input while C P super H iak.es 
ljp^{Ms) as input. Therefore we use the relation 

Mm,) = MM<,) + ^/.log(^). (58) 

The difference between using DR or on-shell quantities as input to the shifts Ap is of 
higher order. 



5.2 Simple approximation of parameter shifts 

It is useful to find a simple approximation for the 0{as) contribution to the parameter shifts 
Ap. It turns out that the shifts in and M^^ in general are less numerically significant 
than the shifts in At. Therefore we investigate an approximate treatment in which the shifts 
in Ml and M^^ are neglected. Furthermore, since we have neglected the D-terms in the stop 
mass matrix and use M^^^y{^s) = ^Y^i^s) = M~^J^{Ms), the stop mixing matrix defined 
in eq. (6) has the simple form 



1 



Accordingly, the relation between |X, 



on— shell | 
t I 



1 

-DR 



and \XY^{Ms)\ simplifies to 



where 



Am? 



Anit 



Am? 

12 



9i 



9i 



on— shell i 



DRi 



Amt \ 
m ) 



2my 



(Am|-Am?J 




6 log 



A'ren 



10 



2mg log 



m^ 



- 1 



fiRe Bo[ml,ml,m 



2 I 2 
m^ + m^ 



and Bo = [-Bo]*^'^) with the scalar integral Bq defined as 

(2vr/i)^-^ 



l)*2mgmf cos(0M3 - 4>Xt), 



(g2 — mg + ie){{q + piY — mf + ie) 



(59) 

(60) 

(61) 
(62) 
(63) 

(64) 
(65) 

(66) 
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Figure 6: Impact of parameter conversions on the neutral Higgs masses M/^, M/^j, M^^. 
The predictions for the Higgs masses, evaluated with the program FeynHiggs, are shown 
as a function of tan/3 in the CPX scenario. Green (dashed): Result where the numerical 
values of the input parameters A^^{Ms), M^^(Ms), M?^(Ms) are directly inserted into 
FeynHiggs without a parameter conversion (i.e., Ap = 0). Blue (solid): Output of FeynHiggs 
if on-shell values of the parameters are inserted that are obtained from the DR input values 
using the full expressions for Ap given in eqs. (48)-(53) at 0{as). Orange (dotted): Output 
of FeynHiggs if the parameter conversion is calculated using the approximation described in 
Sect. 5.2. 



5.3 Numerical examples in the CPX scheme 

In Fig. 6 we investigate the impact of the 0{as) parameter conversions on the predictions 
for the neutral Higgs masses M/^, Mh2, M^^ (see Sect. 10.3 for a discussion of the 0{at) 
terms). We compare the case where the numerical values of the input parameters A^^{Ms)^ 
mW{Ms), MP^{Ms) in the CPX scheme (for Mh± = 140 GeV) are directly inserted as 
input into the program FeynHiggs (in which the parameters are interpreted as on-shell quan- 
tities) with the case where a proper conversion of the DR input values to on-shell parameters 
has been carried out at 0{as) according to eqs. (48)-(53). One can see that the parameter 
shifts have a very large numerical impact on the prediction of the lightest Higgs mass of 
more than 20 GeV in the region of small tan P, while the corresponding effects on the pre- 
dictions for and Mh^ are typically below the GeV level. The result for M/^^ indicates 
the well-known fact that corrections of 0{at<ys) in the MSSM Higgs sector can be numeri- 
cally very important. The numerical effects found here in the CPX scenario are even larger 
than the corresponding shifts in the case of real parameters as discussed in Ref. [77]. As a 
consequence, it is obvious that a proper conversion of parameters at least for the prediction 
of the lightest Higgs mass is indispensable for a meaningful comparison of results obtained 
in different renormalisation schemes. 

Also shown in Fig. 6 is the result obtained from employing the approximate parameter 
conversion as given in Sect. 5.2. One can see that the result for the approximate treatment 
is close to the one obtained with the full parameter conversion. This indicates that the main 
impact of the parameter conversion is indeed caused by the shift in the absolute value of 
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the trilinear coupling AjAtl, as expected from the discussion above. From the expressions 
given in Sect. 5.2 one can furthermore see that there is a significant dependence on the phase 
{(pMa — 4>Xt) s^iid that the gluino mass IM3I plays an important role. We will use the full 
0{as) parameter conversions throughout the rest of the paper. However, we note that this 
approximate parameter conversion can be useful in situations where the inclusion of the full 
0{as) parameter conversions is impractical. 



5.4 Reparametrisation of rrit in the neutral Higgs self-energies 

The difference between parametrising the neutral Higgs self-energies in terms of the on- 
shell top mass and parametrising the neutral Higgs self-energies in terms of the MS top 
mass is formally a three- loop effect. Previously, we have chosen to use an on-shell top 
mass. We shall now investigate the numerical effect of parametrising in terms of the MS top 
mass rut := m^^'^'^ (rrit). In order to simplify the following discussion, we shall assume no 
resummation of tan /3-enhanced terms. 

So far, we have been using neutral Higgs self-energies expressed in terms of the on-shell 
top mass: 

= Syi(rnr) + sllUuk(^r) + S?i(mr) + /^.o.t. (67) 

where '/i.o.t.' stands for 'higher order terms', and the on-shell top mass and the MS top mass 
are related through 

rrit = my + OgX 
4 

X = -—mt (68) 
on 

The difference between using an on-shell top mass and the MS top mass in x is not important, 
since it will affect the calculation only at the 4-loop level. (We shall use mf^). 
We substitute for mf^ in J]{mf^) 

t{mf') = t^^l{mt-asx) + ±^^^_y^^{mt-asx) + ±P^^{^^^ (69) 
and perform an expansion around a; = 0, keeping terms up to 0{atas) 

±imf') = sS,(mO-«sxSS;(m,) + sii_Y,k(m^^ (70) 

The reparametrisation according to eq. (68) has generated an additional term —asxTiy^^^m), 
which is of two-loop order. In order to calculate this additional term, we require the explicit 
expressions for the leading 1-loop Yukawa corrections to the neutral Higgs self-energies, 
which we give explicitly in eq. (97)-eq. (102) below. 

We first substitute for the stop masses their expressions in terms of the top mass and the 
soft SUSY-breaking parameters, in order to ensure that the top mass dependence is explicit 
everywhere. We employ here the Yukawa approximation, i.e. eq. (5) simplifies to 

1 



Ml + MlT\[Ml-MlY + Aml\X,\^ . (71) 
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Figure 7: Neutral Higgs masses in the CPX scenario at tan/3 = 11 and Mh± = 140 GeV 
as function of phase = (f>At = (pMa- Blue (dashed): the Higgs self-energy calculation is 
parameterised in terms of mf^, orange (solid): the calculation is parametrised in terms of 
rut- 

Then we make the substitution nit ^ rat — oisX and expand around x = to obtain terms at 
0{atas). We then edit the FeynHiggs code to include these additional terms and make use 
of the option in FeynHiggs where the tree level stop sector parameters are calculated using 
rrit rather than mf^.^ 

5.4.1 Numerical results 

Fig. 7 shows the neutral Higgs masses in the CPX scenario at tan /3 = 11 and Mh± = 
140 GeV as function of phase = 0At = 0M3, with the Higgs self-energy calculation pa- 
rameterised in terms of m^'^ (blue, dashed) and in terms of (orange, solid). These 
results include the resummation of tan /3-enhanced terms. We can see that, even in the 
CP-conserving limit, parameterising the calculation in terms of Wit rather than mf^ can 
increase the lightest Higgs mass by 5.3 GeV. However, note that we will be interested in 
the maximally CP- violating case = n/2, where we can see that the effect on M/^ is more 
modest: a 1.6 GeV absolute increase or a 2% relative increase. 



6 Higgs cascade decay 

We use the general expression for 2-body neutral Higgs boson decays^, 

TiK^hK) = (72) 

where Ai is the matrix element. 

^Note that the method described here differs from the FeynHiggs option runningMT = 1 at the 2-loop 
order. 

^In this convention, the capital letter F denotes a decay width when it has an argument which explicitly 
contains the symbol (e.g. r{ha — > hbhc)) and a vertex function when it does not (e.g. F^^'v). 
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Figure 8: The leading vertex corrections to the decay hi hjhk, involving t/t loops, (x, y, z 
= 1,2) 



For identical final state particles, hh = 



Mh 4M2 



and the symmetry factor 5 is |. For the case hb 7^ he 



P = ^^J^'^t + K + K - 2 {MlMl + MlMl + MlMl), (74) 

and the symmetry factor is 1. 

Since the lowest order contribution involves only scalar particles, the tree level decay 
width has a very simple form. For example, the h — )■ AA tree level decay width is given by 



6.1 Calculation of the genuine hi — > hjhk vertex contributions 

We calculate the full IPI (one-particle irreducible) 1-loop vertex corrections to the hi — )■ hjhk 
decay width within the Feynman-diagrammatic approach, taking into account the phases of 
all supersymmetric parameters, hi, hj, hk are some combination of the tree level Higgs fields 
h,H,A. The programs FeynArts [119-121] and FormCalc [121,122] are used to draw and 
evaluate the Feynman diagrams using dimensional reduction, and LoopTools [122] is used 
to evaluate the majority of the integrals. We use rrib = m^^'^^ {m^^) and a top mass of 
rrit = 'm^^'^'^'\mf^) = mf'^/ (l + ^^^(mP'^)) in the t, i, b, b masses which enter the genuine 
vertex corrections in order to absorb some of the higher order SM QCD corrections. We use 
a unit CKM matrix and assume no squark generation mixing. 



6.1.1 Leading corrections (Yukawa terms) 

At low to moderate values of tan/3, the leading corrections to the hi — )■ hjhk vertex are the 
Yukawa terms from the t, t sector. These arise from the diagrams shown in Fig. 8. 

As a first step, we select only terms proportional to m^/{M^s^) (Yukawa terms') and 
perform the calculation at zero incoming momentum i.e. T^ijip^ = 0)- this way, we obtain 
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compact analytical expressions for the leading corrections to the hihjhk vertex. We find that 
there are no counterterms contributing to the hihjh^ vertex in this approximation. 

Note that, for consistency, the stop masses m^^ and m^^ must also be calculated in the 
Yukawa approximation according to eq. (71). We therefore arrive at the following expressions 
for the leading Yukawa corrections in the t/t sector, which we can express as corrections to 
an effective coupling iX^^ = iA*''"'^ + iAX^^^. 

For vertices involving the CP-even tree-level Higgs bosons for the case m^^ ^ m^^ (the 
expressions for m^-^ = m^^ are given in App. C): 
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where 



Yt = At + tpfi* 
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and fn'^^^ are the stop masses in the Yukawa approximation, as given by eq. (71). Cii2-i22; 
^^1112-1222 and ^^11122-11222 are functions of Cq, Dq and Eq scalar integrals, respectively. Since 
we are describing a process with 3 external legs, Dq and Eq do not appear explicitly in the 
Feynman diagrams. However, these functions are very useful for simplifying the vertex 
expressions. 

The 1-loop corrections to a hihjhk vertex involving at least one CP-odd eigenstate (again, 
for m^-^ 7^ m^^) are given by 
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These compact, momentum independent expressions have the advantage that they are 
extremely easy to implement into a computer code. In this form, we are also able to see that, 
despite including the effect of complex phases, these corrections are themselves entirely real. 

In order to convert these corrections to the /i, A basis, we use the mixing matrix from 
eq. (2). For example, A^,^^ = X^A - So^^X^^aa + c.AA^^^^. 

In the MSSM with real parameters and m^-^ ^ m^^; these corrections reduce to the form 
(here we drop the subscripts on C112-122, ^^1112-1222, ^^11122-11222 for brevity) 
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and AA^^^;^ = AA^^^;^ = AA^^^;^ = AA^^J^ = 0. These expressions correspond 
exactly to the results for the leading Yukawa corrections to the triple Higgs vertex in the 
MSSM with real parameters published in Ref. [15]. 

For completeness, we also give the leading Yukawa corrections to the neutral Higgs self- 
energies in the complex MSSM, which we used when investigating the effect of reparametris- 
ing the Higgs self-energies in terms of the MS top mass, as described in Sect. 5.4. These 
expressions can be found by considering diagrams involving t, t, b, b loops only and selecting 
those terms proportional to mf/My^^ ('Yukawa terms'). The resulting corrections will be 
finite and proportional to mf. As before, the renormalisation constants Sti3, SM^, (5M| 
and SZij, are all zero in this approximation and the incoming momentum is taken to be 
zero. These expressions also involve the Higgs tadpoles and the charged Higgs self-energy 
(since Mh± is the input parameter rather than Ma), which is also taken at zero incoming 
momentum, such that 

^H-H+ (0) . (96) 

The leading corrections to the renormalised neutral Higgs energies in the Yukawa approxi- 
mation for nii^ 7^ are thus given by [109]: 
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Ci2L = Co(0,0,0,m|,m|,M9. (103) 

Note that the Co integrals do not appear automatically, as no 3-point functions are calculated. 
However, substituting Cq integrals for combinations of the Aq and Bq integrals which appear 
naturally in the calculation (all at zero momentum) does make the self-energy expressions 
more compact. The expressions in the case where m^-^ = m^^ are given in App. C. 



6.1.2 Full 1-loop IPI vertex corrections 

For the full IPI 1-loop corrections to the hihjhk vertex, we need the relevant counterterms. 
Note that, for triple Higgs vertices with an external Higgs boson A, the field renormalisation 
constant 6Zag is required in order to ensure that the vertex is UV-finite. We have extended 
the FeynArts [119-121] model files in order to include these counterterms. 

Examples of Feynman diagrams contributing to these vertex corrections are shown in 
Fig. 9. 

We also investigated the effect of including loop-corrected masses and couplings of the 
Higgs bosons in the one-loop contributions to the hihjh^ vertex, instead of the tree level 
quantities. In order to ensure the UV divergences cancelled, we transformed the couplings of 
the internal Higgs to the other particles using a unitary approximation to the matrix Z (by 
setting the momentum in the neutral Higgs self-energies T^h-h^ to {rn\. + m\.)/2), which we 
implemented into a FeynArts [119-121] model file. For consistency, the loop corrected Higgs 
masses of these internal Higgs bosons were also calculated using this unitary rotation matrix. 
(Note that we continued to use the full Higgs masses and Higgs propagator corrections for 
the external Higgs bosons). These corrections were numerically insignificant in the examples 
investigated. 

6.2 Combining the IPI vertex corrections with propagator cor- 
rections to obtain the full ha — > hbhc decay width 

We can combine vertices involving the tree level Higgs bosons hi,hj,hk with the wave- 
function normalisation factors contained in the matrix Z, which contain self-energies from 
the program FeynHiggs, in order to obtain processes involving the loop-corrected states 
ha, hfj, he as the external particles (as discussed in Sect. 4.2). 
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Figure 9: Examples of generic diagrams (showing only one of the topologies) contributing 
to the processes hi — )■ hjhk. hi, hj, h^ are the physical Higgs fields at tree level {h, H,A), f 
are SM fermions, / are their superparters, X ^i-re neutralinos and charginos, V are vector 
bosons, H denote the neutral and charged Higgs bosons and the Goldstone bosons, u are 
Faddeev-Popov ghost fields. 



These 'Z-factors' can be used in conjunction with tree level^ vertices ^^^^h^h^ using (sum 
over k) 

The decay width is then given by eq. (72) with \M.\^ = l-^l^ = \K\hf- 

Note that this means that our decay width will contain pieces of type (1-loop) x (1-loop), 
which is necessary since the tree level coupling is often smaller than the leading loop correc- 
tions to the coupling. 

We obtain our full result by combining the complete genuine 1-loop vertex corrections 
and the corrections involving 1-loop Goldstone and Z boson self-energy contribu- 
te, Zse 



tions ^iT^^hjhk "with the Z-factors, such that (sum over k) 

■ptrec , -plPI.l-loop / ,^2 7i^2 /\/f 2 \ , "pG.Z se / 



"pfull ry ^ ^ 



(105) 

The genuine 1-loop vertex corrections ^^^1%^°°^ contain the full momentum dependence 
and therefore depend on the loop-corrected masses M^^, at the external legs. How- 

ever, as discussed in detail in Sect. 4.3, unphysical poles from diagrams involving Goldstone 
and Z boson self-energies can be avoided by approximating the external momenta to the tree 



level values in the corresponding contributions, i.e. F^,'^^*^^^ is a function of 
Again, the decay width is then given by eq. (72) but with = |F^"jJ,^^^p. 



^This method can also be used for the effective vertices ^hlhbha ^-^^ ^lilh^ha ~ ^ck^bj^ai^h^hjh^ 
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6.3 Numerical Results 



6.3.1 /i2 — > hihi decay width 

We will now investigate the importance of the full 1-loop genuine corrections through their 
numerical impact on the h2 — ?■ hihi decay width. All the results plotted in this section 
include the wave-function normalisation factors, through the matrix Z. The case where 
only wave-function normalisation factors but no genuine one-loop vertex contributions are 
included will be denoted 'tree'. 

Fig. 10 compares the 'tree' result with the full result which includes the genuine vertex 
correction and all propagator corrections, as described by eq. (105), in the CPX scenario. In 
Fig. 10 (left), we can see that 'tree' and full decay widths are very different. The full result 
has a peak (i.e. local maximum) at tan/3 = 8.7. There is a corresponding peak in the 'tree' 
result at tan /3 = 5.5. However this peak is 7.5 times smaller than the peak in the full result. 
The sharp increase in the full result at low tan/3 is because we have chosen to keep M^^ 
constant, which requires a rapid change in Mh± in this part of parameter space (the 'tree' 
result also exhibits this behaviour, but at tan/3 ~ 2.5, which is not shown on the graph). 

In Fig. 10 (right), we can see that both the 'tree' and full decay widths decrease as the 
lightest Higgs mass M/^^ increases, as we would expect from the kinematics. Again, the 
'tree' level result is heavily suppressed compared to the full result. We can conclude that 
calculations of triple-Higgs couplings which just combine the propagator corrections with the 
tree level vertex but do not take into account genuine vertex corrections are an extremely 
poor approximation to the full result. 

Fig. 11 demonstrates the pronounced dependence of the results on the complex phase 0, 
where = (pAt = 4>Ab = (pA^ = (pMs, at tan/3 = 11, Mh± = 300 GeV (all other parameters 
are taken from the CPX scenario). We can see once again that while the 'tree' result gives 
qualitatively similar behaviour, the peaks are less than a quarter of the peaks in the full 
result, and the positions of the troughs are shifted. It is interesting to note that the genuine 
vertex corrections play a highly significant role over essentially the full range of 0, including 
the points at = and = tt where the parameters are entirely real. 

Fig. 10 and Fig. 11 also demonstrate the results of two methods for obtaining effective 
hihihi couplings: the leading Yukawa corrections eqs. (76)-(89) and the fermion/sfermion 
corrections taken at = 0. Both approximations are a big improvement over the 'tree' 
result. 

In Fig. 10, at the peak at tan /3 = 8.7, the result using /, / at = is within 14% of 
the full result and the 'leading Yukawa' approximation is within 27% of the full result. As 
we can see from Fig. 11 at = 0, in some parts of parameter space, the 'leading Yukawa' 
approximation accidentally performs better than the /, / at = terms, although it is a 
less complete approximation to the full result. 

This calculation of the /i2 — > hihi decay width also applies to the MSSM with real 
parameters. In Fig. 12, we show the results for the decay width r{H — )■ hh) in the m™^^ 
scenario as a function of tan /3 at Mh± = 300 GeV. As we saw in the CPX scenario, 
using the tree level triple Higgs vertices combined with the propagator corrections is a 
poor approximation to the full genuine triple Higgs vertex corrections combined with the 
propagator corrections. The 'leading Yukawa' approximation and using the fermion/ sfermion 
corrections at = in the genuine vertex corrections are far better approximations to the 
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Figure 10: The decay width r{h2 — hihi) in the CPX scenario as a function of tan/3 at 
M/ij = 40 GeV (left) and as a function of M/^ at t/3 = 6 (right). Mh± is adjusted to give the 
M/jj as required. All results include the propagator corrections, 'tree' indicates that the tree 
level triple Higgs vertex has been used, 'leading Yukawa' includes the leading genuine vertex 
corrections given in eqs. (76)-(89). 'f,sf p^=0' includes contributions to the genuine vertex 
corrections from fermions and sfermions only and approximates the external momenta to 
zero. The curve labelled 'full' shows the result including the full genuine vertex corrections. 



full result. 

6.3.2 Effective coupling approximation for the lightest neutral Higgs boson 

As we have seen, the effective triple Higgs vertices obtained using the leading corrections 
in the Yukawa approximation (as given in eqs. (76)-(89)) and the effective triple Higgs 
vertices obtained using the fermion/sfermion corrections at = in the genuine vertex 
corrections and combining with the propagator corrections have both performed rather well 
as approximations to the /i2 — > /ii/ii or if — )■ hh decay width. We have also seen that 
it is possible to get a large enhancement in these decay widths from the genuine vertex 
corrections. We shall now apply these approximations to the triple coupling of the lightest 
Higgs boson. 

Fig. 13 shows this effective coupling squared, supplemented by propagator corrections 
and normalised to the tree level SM value (with equal Higgs mass). We can see that there 
is a suppression with respect to the tree level SM value if no genuine vertex corrections are 
included. This holds even in the limit of large rriA- This is because the SM tree-level coupling 
involves the square of the physical Higgs mass whereas the effective coupling in this limit 
involves the square of the tree-level mass of the lightest MSSM Higgs boson. 

In Fig. 13, we also show the areas which are already excluded at the 95% CL by the 
LEP Higgs searches, which we have determined using the method described in Sect. 10.2 
below. Including the genuine vertex corrections gives an overall enhancement of the effective 
couplings squared of approximately 1.2 to 1.6 in the parts of this parameter space which are 
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Figure 11: The dependence of the h2 — )■ hihi decay width on the phase (f) = (pAt = (pA^ = 
(j)A^ = 4>M-i (CPX scenario at tan/3 = 11, Mh± = 300 GeV). The full result is compared 
with different approximations, which are the same as those specified in Fig. 10. 

not yet excluded. This could have interesting implications for the sensitivity of searches at 
the LHC and LC to effects of the triple-Higgs coupling in the MSSM. 



7 Higgs decay to SM fermions 

7.1 Calculation of the ha — > bb decay width 

We use the general expression for two-body decay widths given in eq. (72) to obtain: 



where N^. is the number of colours. The mass dependence of the squared matrix element 
will be affected by the CP properties of the Higgs boson. For example, at lowest order. 



T,MK^m - ^^(l-gVlA-r. (107) 



where \X]^;;^\^ = Xf'^'^sl, \X%2\^ = X^^^'^cl and \X'lll\^ = A[°^'^4 with x = 3/2 for the CP 



even states h, H and x = 1/2 for the CP odd state A. X^^^ was defined in eq. (37) 
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Figure 12: The decay width r{H — hh) in the m^^^ scenario as a function of tan /3 at 
Mh± = 300 GeV. The full result is compared with different approximations, which are the 
same as those specified in Fig. 10. 



7.1.1 Standard Model QED corrections 

The real and virtual QED contributions to the Standard Model — )■ bb decay width lead 
to the 1-loop correction 
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for M^^'^ ^ rrif, as derived in Ref. [123]. In this limit, this equation can be used for the 
QED corrections for both the scalar and pseudoscalar MSSM Higgs bosons [31]. We will 
therefore use the correction term 
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(109) 



in our MSSM calculation. 



7.1.2 Standard Model QCD corrections 

If the factor Qla in eq. (108) is replaced by the factor C/a,(M|^), the expression for the 
1-loop QCD correction to the — )■ bb decay in the Standard Model is obtained, as shown 
in Ref. [123]. Cf is the quadratic Casimir operator (C/ = |) and ^^(M^^^) is the running 
coupling. Including the tree level result gives 
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Figure 13: The triple coupling of the lightest Higgs boson squared, normalised to the SM 
tree-level value (with equal Higgs mass). Left: m™^^ scenario at Mh± = 300 GeV, varying 
tan/3. Right: CPX scenario, tp = 11, Mh± = 300 GeV, varying the phase (p = (pAt = 
= 4>Ar = 'pM-i- AH results include the propagator corrections, 'tree' indicates that the 
tree level triple Higgs vertex has been used, 'leading Yukawa' includes the leading genuine 
vertex corrections using the Yukawa approximation, as given in eqs. (76)-(89). 'f,sf p^=0' 
includes contributions to the genuine vertex corrections from fermions and sfermions only 
and approximates the external momenta to zero. Yellow (shaded): excluded at 95% CL by 
the LEP Higgs searches (see Sect. 10.2 for details). 



(110) 

where we have factored out the Yukawa coupling from the term in the square bracket. 

In the mass range we are interested in, Mf^ ^ rUb, this leading logarithmic approxi- 
mation is not sufficient. However, the large logarithmic corrections can be absorbed into a 
running bottom quark mass. Substituting the relation between the on-shell b-quark mass 
and the running SM MS b-quark mass into eq. (110) gives 
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and we can choose /iren = Mff in order to cancel the logarithmic term. 

In practice, we parametrise our calculation in terms of m^^'^^^Mh^). Therefore, in order 
to encompass the full 1-loop Standard Model-type QCD corrections in our calculation, we 
will need to add a correction 
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to the ha — bb decay width. 
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Our method differs from that of Ref. [30], which includes some higher order terms in 

asifJ'ren) "^fo^^''^*^ (/^ren) and an extra term proportional to as{ml) / as{MlJ . Our method 
also differs from Ref. [94], which includes terms proportional to a^. However, some of these 
terms depend on the CP properties of the Higgs, and thus are not trivially extendable to 
the complex MSSM. Both Ref. [30] and Ref. [94] restrict their analyses to the MSSM with 
real parameters. 

7.1.3 Full 1-loop IPI hi — 7- bb vertex corrections 

In order to calculate the 1-particle irreducible vertex corrections at 1-loop, we have once again 
extended the FeynArts model file to include the relevant counterterms. We have evaluated 
the complete contributions in the MSSM, apart from the QED- and QCD-type corrections, 
for which we use the results given in eq. (109) and eq. (112), respectively. (Note that we 
include the photon contribution to the W boson mass counterterm) . We include all complex 

phases. As discussed above, we use rri^^'^^\Mh^) in these corrections in order to absorb 
some of the higher order terms. As before, we use a unit CKM matrix and include no squark 
generation mixing. 

7.1.4 Resummed Amf, corrections to hi bb 

In order to resum the leading SUSY QCD (and higgsino) corrections for the case of large 
tau/S in the limit of heavy SUSY particles, we use the effective vertices v^^^ from the 
effective Langrangian in eq. (33). 
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When combining this contribution with the full genuine vertex corrections, we need to 
avoid double-counting of the 1-loop corrections involving gluinos or higgsinos. Therefore, we 
subtract the 1-loop part of these effective vertices, i.e. we use effective couplings of the form 
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7.1.5 Combining these contributions with propagator corrections to obtain the 
full ha — )■ bb decay width 

The tree level, 1-Ioop IPI hi — j- bb vertex function and the additional Am^ corrections are 
then combined with the propagator corrections from the hihj self-energies (via Z) and from 
the hiG and hiZ self-energies: 



^habb 



The arguments to r!;^,^^ and T^''^t^'^ denote the external momenta used. F^™;'"'^"^ 

hibb hibb hibb 



loop 



denotes the use of (^f^f^'' — ^Vsb '^ ^°°^) " ^habl is combined with the external fermion wave- 
functions, then we take the squared modulus and sum over external spins to get \M.habl\^- 
The full ha — )■ bb decay width is thus found using 



T'-\K^bb) = [l + '^QCD + 5QED]^^^yi-^|-M.>.|% (118) 

which is an extension of the method used to combine QED, QCD and Z-factor corrections 
in Ref. [30]. 



7.1.6 Numerical Results 

Fig. 14 illustrates the decay widths T{hi — )■ bb) (left), F(/i2 — >■ bb) (centre) and F(/i3 — )• bb) 
(right) as a function of the charged Higgs mass for the CPX scenario with tan /3 = 20. All 
results include the propagator corrections, incorporated via the matrix Z, and the Goldstone 
and gauge boson mixing contribution. We note that the hi and /12 decay widths have steep 
gradients at Mh± ~ 157 GeV due to a 'cross-over' effect in the masses (i.e. M^^ and M^^ 
approach each other). At Mfj± ~ 150 GeV, hi is mostly A, /12 is mostly h and /13 is mostly 
H, whereas at Mfj± ~ 200 GeV, hi is mostly h, h2 is mostly A and h^ is mostly H. 

The impact of including the resummed Ami, contributions is very significant, as we 
can see from comparison of the result using the propagator, SM QCD and QED corrections 
CQED, QCD') and the result when the resummed Amj, contributions are included in addition 
('QED, QCD, Amb)- Fig. 14 also includes the full decay widths ('full'), which combine the 
propagator, QED, SM QCD, SUSY QCD corrections with the full electroweak genuine vertex 
corrections, as described by eq. (118), and the effect of including the full corrections turns 
out to be relatively small but not negligible. At Mh± = 600 GeV, these extra electroweak 
corrections are 5.6%, 5.7% and 5.0% for the decay of hi, /12, hs, respectively. 

It is possible to approximate the 'QED, QCD, Am^' result using an effective 6-quark 
mass of m^^'^ ^^rrit) / \1 + Arribl (with no additional Arub contributions), which we call 'QED, 
QCD, effective nib - The resulting decay widths are also shown in Fig. 14. In general, this 
approximation performs well. However, note that in the decoupling limit, where the hi is 
SM-like, the effect of the SUSY QCD corrections should cancel out in the hi decay width. In 
the approximate 'QED, QCD, effective rrib result, this does not occur, giving a offset which is 
significant in relative terms, although small in absolute terms. In the region Mh± = 200 GeV 
to Mh± = 600 GeV shown here, this offset between the 'QED, QCD, effective nib result 
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Figure 14: The decay widths of the neutral Higgs bosons to two b-quarks in the CPX scenario 
with tan f3 = 20. See text for explanation of the various approximations used. 

and the approximate result is more than 44% of the 'QED, QCD, Arrib result, while the 
absolute difference is less than 0.002 GeV. 

The ha — )■ bb decay width in the CPX scenario is highly dependent on (pM-j- Fig. 15 shows 
the decay widths of the two lightest Higgs bosons into b-quarks at (pMa = 0, (j)M3 = 7r/2 and 
0M3 = 37r/4. Whereas the inclusion of the resummed Am^-type corrections suppresses the 
decay width at (pMs = 0, at (pM^ = 7r/2 and (pMs = 37r/4 these corrections have caused an 
enhancement. The value of Arrib in these plots is 1.05 — 0.12^, —1.162,-0.75 — 0.86i for 
=0, 4>M3 =7^/2 and (pM^ = 37r/4, respectively. 

7.2 Calculation of the ha — > r+r decay width 

The calculation of the ha — t- r+r^ decay width is similar to that of the ha — )■ bb decay 
width, with the simplification that no QCD corrections are required. We have calculated the 
full 1-loop genuine vertex corrections and supplemented these with propagator corrections 
(including 1-loop mixing with Goldstone and Z bosons) and QED corrections. As before, we 
have included all complex phases. 

8 Higgs cascade decay branching ratios 

Accurate predictions for Higgs branching ratios are vital for Higgs phenomenology. In partic- 
ular, they are required as part of calculations of cross sections of collider processes involving 
the production and decay of an on-shell Higgs boson, which are often performed using the 
narrow width approximation. In Sect. 10.2, we will use Higgs branching ratios for the CPX 
scenario in conjunction with the LEP topological cross section limits. In order to understand 
the resulting exclusions, it will be necessary to refer to the behaviour of the contributing 
branching ratios. 

We combine the ha — )■ hbhc decay widths calculated in Sect. 6 with the ha — )■ bb and ha — )■ 
T~T~^ decay widths calculated in Sect. 7. As we have discussed, these decay widths include 
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Figure 15: The decay widths of the two hghtest neutral Higgs bosons to two b-quarks in the 
CPX scenario with tan/3 = 20 at various values of (/)m3- Solid: ha = hi, dashed: ha = h2- 

the full 1-loop genuine vertex corrections and are combined with propagator corrections^ 
obtained using neutral Higgs self-energies from the program FeynHiggs [6-10], which include 
the leading 2-loop contributions. The 1-loop propagator mixing with Goldstone and Z bosons 
is also consistently incorporated. These results take into account the full phase dependence 
of the supersymmetric parameters. For the ha — )■ bb decay width, the Am^ corrections are 
resummed in a way that preserves the phase dependence. We take all other decay widths 
from the program FeynHiggs [6-10] (these decay modes are subdominant in most of the 
regions of MSSM parameter space). 

Fig. 16 (left plot) illustrates the pronounced dependence of the /i2 — )■ hihi branching 
ratio on tan/3 and Mh^. We see that this decay mode is significant and often dominant in 
almost all of the regions where it is kinematically allowed. We can see that the characteristics 
of the /i2 — hihi branching ratio are largely determined by the behaviour of the /i2 hihi 
decay widths (in Fig. 10 we showed this decay width for two slices of CPX parameter space). 
Note, in particular, the narrow 'knife-edge' region of very low /12 hihi branching ratio, 
which occurs at tan /3 ~ 4.5, where the /i2 — ?■ hihi decay width tends to zero. The behaviour 
of the h2 — )■ bb branching ratio is also heavily dependent on the h2 — )■ hihi decay width 
where the /i2 — hihi decay is allowed kinematically, since in this region, the /i2 hihi 
decay usually makes up the majority of the total decay width. Over the majority of CPX 
parameter space, BR(/i2 hihi)^BR{h2 66)-hBR(/i2 r+r~)~ 1, and BR{h2 -> r+r") 
is comparatively small. 

^Notice also that there are some points within the CPX parameter space that are shown here without a 
branching ratio value and that the edge of the allowed parameter region is uneven. These are points where 
either the mass calculation or the Z-factor calculation did not produce a stable result because the terms 
involving double derivatives of self-energies were non- negligible, as described in eq. (17). 
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Figure 16: The branching ratios of neutral Higgs bosons into other neutral Higgs bosons in 
the CPX scenario. 

The Higgs cascade decays for the heaviest neutral Higgs also dominate in the majority 
of the region where they are kinematically allowed. The — )■ hihi branching ratio (middle 
plot of Fig. 16) also has a narrow region at tan/3 ~ 4 — 5 in which, while the — )■ hihi 
decay is kinematically allowed, the decay width is nevertheless suppressed, characteristically 
similar to the suppressed region we observed in the /12 hihi branching ratio. In particular, 
topologies involving /13 can be relevant to the LEP exclusions in the region 10 ^ tan/3 ^ 30, 
^ 60 GeV for variations of the CPX scenario. In this region of parameter space, the 
hs — hihi decay width is a crucial contribution also to the — >■ bb branching ratio. The 
^3 ~^ h2hi decay width (the /13 — >■ /i2^i branching ratio is shown in the right plot of Fig. 16), 
on the other hand, is of less phenomenological interest in this scenario, since this decay width 
dominates the total width in a region which, as we will see, can be excluded at the 95% 
CL using limits on the cross sections of topologies involving hi and /i2- 
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9 Normalised e^e Higgsstrahlung and pair production 
cross sections 

In order to examine the effect of our predictions for the Higgs branching ratios on the size of 
the regions of CPX parameter space which can be excluded by the LEP Higgs searches, we 
need to consider the Higgsstrahlung and Higgs pair production cross sections, normalised to 
a reference cross section cXref • In the MSSM with real parameters, the corrections to the LEP 
Higgsstrahlung and LEP pair production processes have been studied in detail [34-39]. In the 
MSSM with CP violation, there has also been considerable interest in accurate predictions 
of these cross sections, although the full 1-loop corrections are not yet available [40-44]. 

For the Higgsstrahlung topologies e^e~ — haZ, the reference cross section a^ci is the 
tree level SM cross section for the process e^e~ — > Z — )■ HZ, for a SM-like Higgs of mass 
Mh = Mh^. For the pair production topologies e+e" — )• hahb, the reference cross section a^ci 
is the tree level MSSM cross section for the process e"'"e^ — )■ Z — ?■ h^A^, where the MSSM 
coupling factor cos^ (/3 — a) has been divided out and the masses of h'^ and taken as Mh^ 
and Mfi^ respectively. This reference cross section can also be expressed in terms of the 
Standard Model Higgsstrahlung production cross section, 

cTrcf = AcT^^(e+e- ^ HZ), (119) 

where A is a kinematic factor which takes into account the different kinematic dependences 
of the SM Higgsstrahlung and the pair production process, i.e. 

A = \Zj{l2Mys + \zH)/\TH, (120) 
A,, = [l-{M, + Myf/s][l-iM^-Myf/s], (121) 

and if is a SM-like Higgs with mass M^- 

For the majority of our scans, we will calculate these normalised cross sections using 
an effective coupling for the ha-Z-Z or ha-h^-Z vertex, which incorporates external Higgs 
propagator corrections. However, we will also examine the effect of including the complete 
t, t, 6, h 1-loop corrections (involving also genuine vertex corrections) in these cross sections. 

9.1 Normalised effective Higgs couplings to gauge bosons 

The matrix Z can be used to create a normalised effective coupling between neutral Higgs 
bosons and Z bosons which takes the corrections to the external Higgs propagators into 
account, though the relations 

Q'haZZ = '^aiOhiZZi (122) 
dh^HZ = '^bj'Z^aigt'^hjZ-: (123) 

where g)l^zz ^^^^ normalised to the SM coupling, such that {g^z'zY — ^in^ (/? — a), (qh'zz)'^ = 
cos^ (/3 — a) and (Qa'zz)'^ — 0- addition, we normalise the g^^^jZ ^^'^^ that {OhAz)'^ — 
cos^ (/3 — a) and (Qh^Xz)'^ — ^in^ (/3 — a). All other g^l^^.z are zero. 

Fig. 17 illustrates the normalised squared effective Higgs couplings to gauge bosons 
Igf^zzl"^ in the CPX scenario. We can see that the hi-Z-Z coupling dominates around 
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Figure 17: The effective couplings of the neutral Higgs bosons to two Z bosons, 
l5'/ifzzP;lfl'ft,2ZzP IS'hsZzP' which include the Higgs propagator corrections calculated 
using the matrix Z. 

the edge of the available parameter space, the h^-Z-Z coupling dominates in a region 
Mh^ < 60 GeV and 7 ^ tan/3, and the h2-Z-Z couphng dominates the region in between, 
such that \gf^zz? + Wh^zz? + \9f,zz? ~ 1- Fig. 18 illustrates the behaviour of \gf^h,z?, 
which can be described as \gh'^«zz\'^ ~ IS'/if/i^zP' where ha-, hb, he are all different. (If a unitary 
approximation to the Z matrix is used, as in the LEP Higgs Working Group analysis of the 
CPX scenario [2], these relations become equalities.) 

The normalised Higgsstrahlung and pair production cross sections can then be approxi- 
mated by these effective couplings, i.e 



(T(e+e -^haZ)/a,c{ 
cr(e+e~ hahb)/arei 



\yhaZZ\ 5 

lyhah^zi ■ 



(124) 
(125) 



9.2 Loop corrections to the Higgsstrahlung and pair production 
cross sections 

We now consider 1-loop corrections to the LEP Higgsstrahlung and pair production processes 
involving the particles t, t, b, b. The structure of these diagrams is shown in Fig. 19, where 
loops involving t, t, b, b are indicated by grey circles. In general, lines labelled with X can be 
h, H, A, G, Z, where the resulting diagram conserves CP at the tree level vertices involving 
the gauge and Higgs bosons. We neglect the electron mass and therefore will not include 
any diagrams in which a Higgs boson couples directly to the electrons. We then include the 
Higgs propagator factors for the Higgs on external legs through the Z matrix, as previously. 

The result of including these corrections on the normalised Higgsstrahlung cross section in 
the CPX scenario at M^^ = 40 GeV is given in Fig. 20 (left). The dashed lines show the result 
for the tree level hi-Z-Z vertex combined with Higgs propagator factors, using eq. (122). The 
solid lines show the result when the full t, t, b, b 1-loop corrections are also included. Similarly, 
the dashed lines in Fig. 20 (right) show the normalised pair production cross sections using 
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Figure 18: The effective couplings of two neutral Higgs bosons to a Z boson, \gf^h2z\'^ ^Idf^hazl"^ 
and 1(7/1^/1321^, which include the Higgs propagator corrections calculated using the matrix Z. 
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Figure 19: Additional corrections to the Higgsstrahlung and pair production processes. The 
grey circles indicate loops involving t,t,b,b and Xj,Xj = {h, H, A,G, Z). Note that not 
all of the combinations are physical, since the tree level vertices must be CP-conserving. 
In addition, we neglect the electron mass, therefore no contributions where a Higgs boson 
couples to the electrons are included. 



eq. (123), and the solid lines show the result if the t,t,b,b 1-loop corrections are also used. 
The effect of including the t, t, b, b corrections turns out to be negligible for the Higgsstrahlung 
process. However, for the pair production process, a more sizable effect is visible, leading to 
an increase of the normalised cross section. This is due to the additional Yukawa coupling 
in the genuine vertex corrections of the Higgs pair production process as compared to the 
Higgsstrahlung process, so that a larger enhancement factor is possible in this case. As will 
be discussed in the follwing section, the impact of this kind of corrections on the coverage 
of the LEP Higgs searches in the CPX scenario is nevertheless rather small. 
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Figure 20: The normalised Higgsstrahlung cross sections (left) and normalised pair pro- 
duction cross sections (right) in the CPX scenario at M^^ = 40 GeV. Dashed lines: using 
eq. (122) or eq. (123). Solid lines: as for the dashed lines except that the t,t,b,b 1-loop 
corrections are also included (neglecting the electron mass). 

10 Confronting the Higgs sector predictions with lim- 
its from the LEP Higgs searches 

10.1 LHWG limits on the parameter space of the CPX scenario 

After the LEP programme finished in 2000, the final results from the four LEP collaborations 
(ALEPH [1,124,125], DELPHI [126,127], L3 [128] and OPAL [129,130]) were combined 
and examined for consistency with a background hypothesis and a signal plus background 
hypothesis in a coordinated effort between the LEP Higgs Working Group for Higgs Searches 
and the LEP collaborations (LHWG). The results showed no significant excess of events 
which would indicate the production of a Higgs boson. In the Standard Model, a lower 
bound on the Higgs mass of 114.4 GeV at the 95% confidence level was established [1], while 
restrictions were placed on the available parameter space of a variety of MSSM benchmark 
scenarios [2], including the CPX scenario [5]''. 

For the purposes of the LHWG analysis, two different programs were used to calculate 
Higgs masses and branching ratios in the complex MSSM: FeynHiggs version 2.0 [10] and 
CPH [5], which was a predecessor of the program CPsuperH [11, 12]. These two codes had 
significant differences in the incorporated higher order corrections, and it was necessary to 
perform a conversion between the two sets of input parameters, due to the different renormal- 

^Note that the definition of the CPX scenario used in Refs. [2,5] difi'ers slightly from the definition used 
in the present paper, as discussed in Sect. 3. 
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isation schemes used in the two codes. As explained below, the parameter conversion used 
was an approximation based on a calculation performed in the MSSM with real parameters. 

Separate analyses were performed using FeynHiggs and CPH. In order to combine these 
results, a conservative method was adopted, in which a point in parameter space was regarded 
as excluded only if it was excluded by both the analysis using results from FeynHiggs and the 
analysis using results from CPH. The LHWG analysis resulted in three unexcluded regions 
of CPX parameter space at 95 % CL: 

• (A) 60 GeV < Mh^ and 3 < tan/3 

• (B) 30 GeV < Mh^ < 50 GeV and 3 < tan/3 < 10 

• (C) GeV < Mh, < 10 GeV and 3 < tan/3 < 20 

The results from the separate FeynHiggs and CPH analyses showed substantial differ- 
ences. In particular, the FeynHiggs analysis had a larger unexcluded region of type B, and 
the CPi/ analysis had a larger unexcluded region of type A, while both results showed similar 
unexcluded regions of type C. We shall concentrate on unexcluded regions of type A and 
B in this paper, since constraints other than those from Higgs searches play a role in the 
unexcluded region C (see, e.g., Ref. [131] for a discussion of region C). 

There was an additional complication, since FeynHiggs as yet does not have a reliable 
calculation for the loop corrections to the triple Higgs couplings in the CP-violating MSSM. 
For the purposes of the ^ FeynHiggs^ analysis, the triple Higgs coupling was therefore obtained 
from CPH., and then combined with Higgs masses and other Higgs sector quantities as 
calculated by FeynHiggs [2]. As we will see, higher order corrections to the triple Higgs 
coupling have a great influence on the size, shape and position of the unexcluded region B. 

10.2 The effect of the new Higgs sector corrections on the Umits 
on the parameter space of the CPX scenario 

The LEP Higgs Working Group for Higgs Searches and the LEP collaborations also published 
their combined results in the form of topological cross section limits at 95% CL, which can 
be applied to a wide range of theoretical models. In each of these topologies, the Higgs is 
produced either through Higgsstrahlung or pair production and decays either to b-quarks, 
tau-leptons or via the Higgs cascade decay. To a very good approximation, the kinematic 
distributions of these processes are independent of the CP properties of the Higgs bosons 
involved, as discussed in Ref. [2]. Therefore, the same topological bounds can be used for 
CP-even, CP-odd or mixed CP Higgs bosons. 

In this section, we will use the topological cross section limits from LEP in conjunction 
with updated predictions for the Higgs masses, couplings and branching ratios. In particular, 
we shall be using our full 1-loop diagrammatic calculation for the hi — hjhk decay processes 
with full phase dependence as described in Sect. 6, combined with renormalised neutral Higgs 
self-energies obtained from the current version of FeynHiggs (which includes corrections at 
Ogatas) with full phase dependence). 

In order to utilise the cross section limits, we use the program HiggsBounds [57,58]. As 
input, it requires the Higgs masses, the normalised production cross sections and the branch- 
ing ratios BR(/ia — )■ hbhb), BR(/ia — )■ bb) and BR(/ia — t- r+r~) for each parameter point. We 
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obtain the Higgs masses as described in Sect. 4 and we calculate the Higgs branching ratios 
as described in Sect. 6. We will begin by using LEP Higgsstrahlung production cross sec- 
tions which include both the full propagator corrections and additional t, i, b, b corrections, 
as discussed in Sect. 9.2. Unless otherwise stated, we will investigate the CPX scenario, as 
defined in Sect. 3 and, unless explicitly stated, we do not consider the possible impact of 
theoretical uncertainties from unknown higher order corrections on the exclusion bounds in 
the parameter space. 

HiggsBounds uses the provided Higgs sector predictions to determine which process has 
the highest statistical sensitivity for setting an exclusion limit for each parameter point, 
using the median expected limits based on Monte Carlo simulations with no signal. It then 
compares the theoretical cross section for this particular process with the experimentally 
observed limit for this process. In this way, only one topological limit is used for each 
parameter point, thus ensuring that any resulting exclusion is valid at the 95% CL. 

However, it should be noted that, in general, the dedicated analyses carried out in Ref. [2] 
for specific MSSM benchmark scenarios have a higher exclusion power than analyses using 
HiggsBounds, since, in a dedicated analysis, the information from different search channels 
can be combined. This can, in particular, lead to an improved result in regions of parameter 
space where several channels have similar statistical sensitivities. 

Fig. 21 (left) indicates which channel has the highest sensitivity and therefore which 
channel will be used for each point in CPX parameter space, to determine whether or not 
it is excluded at the 95 % CL. The channel hiZ — )■ bbZ (■) has the highest statistical 
sensitivity at the edge of the CPX parameter space where tan /3 is low or M/j^ is high, due to 
the fact that the coupling of the lightest Higgs boson to two Z bosons is unsuppressed in this 
region, as we saw in Fig. 17. Similarly, in a band adjacent to this, where the coupling of the 
second heaviest Higgs boson to two Z bosons is unsuppressed, the Higgstrahlung processes 
—J- bbZ (M) and h2Z — > hihiZ — > bbbbZ {O) have the highest statistical sensitivity. 
Finally, in the upper left region of the plot, where Igf^h^zl"^ high, the pair production 
processes /i2^i bbbb{M) and h2hi — )■ hihihi — t- bbbbbb{ ) have the highest statistical 
sensitivity. The part of the parameter space in which the processes directly involving the 
/i2 — )■ hihi decay {h2Z — )■ hihiZ — )■ bbbbZ (D) and h2hi — )■ hihihi — > bbbbbb{ )) dominate 
occurs in a region with an increased /i2 — ^ hihi branching ratio at tan /3 ~ 5 — 10, as we 
saw in Fig. 16 (left), and which is influenced by the peak in the h2 — )■ hihi decay width as 
shown in Fig. 10 (left) at tan /3 ~ 8. Fig. 21 (right) differs from Fig. 21 (left) in that only the 
propagator corrections have been used when calculating the predictions for the LEP Higgs 
production cross sections (i.e. we use the normalised squared effective couplings \g'}^^zz\'^ ^"^^ 
IdKhtzl"^ described in Sect. 9.2). We see that the graphs are very similar, with the main 
difference being the reduced size of the hiZ — t- bbZ (■) region at tan /3 ~ 5, M/j^ ~ 25 GeV. 

In Fig. 22 (left), we have compared our theoretical cross section predictions for each 
parameter point in the CPX scenario with the observed topological cross section limits 
obtained at LEP for the channel with the highest statistical sensitivity at that point, in order 
to obtain exclusions at 95% CL. The Higgsstrahlung topologies where the Higgs decays to 
b-quarks are unable to exclude Higgs masses above Mh-^ ^ 114.4 GeV, as we would expect, 
since this is the limit on the mass of a SM-like Higgs boson [1]. The upper edge of the 
excluded area in the h2Z — bbZ (M) region has a similar shape to the Mh^ ^ 114.4 GeV 
contour and, at tan /3 = 20, occurs at a position AMh^ ~ 4 GeV relative to it. As before, we 
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Figure 21: The coverage of the LEP Higgs searches in the M/^^-tan/? plane of the CPX 
scenario, showing the channels that are predicted to have the highest statistical sensitivity 
for setting an exclusion limit. The colour codings are: red = hiZ — )■ 66Z(B), blue = 
— )■ hhZ{M), white = — )■ hihiZ — ?■ bbbbZ {O), cyan = h2hi — bbbb{M), yellow = 
^2^1 — ^ hihihi — >■ bbbbbb{ ), green = h^hi — )■ bbbb{M), purple = other channels (■). Left: 
full result (i.e. including the extra corrections to the Higgstrahlung and pair production cross 
sections described in Sect. 9.2). Right: using effective coupling approximation for production 
cross sections. 

call the unexcluded area in the top right region of the plot, 'unexcluded region A'. It has a 
narrow 'tail', which extends to lower tan/3, one side of which is bounded by the limit for a 
SM-like M/i2 and one side of which is bounded by the edge of the region where the channel 
hiZ — bbZ (M) has the highest statistical significance, as shown in Fig. 21 (left). Fig. 22 
(right) differs from Fig. 22 (left) in that only the propagator corrections have been used 
when calculating the predictions for the LEP Higgs production cross sections. Once again, 
we see that the extra t, t, b, b corrections only have a very small numerical effect. Therefore, 
we will neglect these corrections in the results that we show below. 



44 



exclusions 



exclusions 




Figure 22: The coverage of the LEP Higgs searches in the M/^^-tan /3 plane of the CPX 
scenario, showing the parameter regions excluded at the 95% C.L. by the topological cross 
section limits obtained at LEP. The colour codings are: green (darker grey) = LEP ex- 
cluded, white = LEP allowed. Left: full result (i.e. including the extra corrections to the 
Higgstrahlung and pair production cross sections described in Sect. 9.2). Right: using effec- 
tive coupling approximation for production cross sections. 
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The plots in Fig. 22 also show an unexcluded region of type B at My,^ ~ 45 GeV and 
tan /3 ~ 8. From comparison with Fig. 16 (left), we can see that the entire unexcluded region 
B in Fig. 22 lies in an area where the h2 — )■ hihi branching ratio is sizable. We examine this 
unexcluded region B in more detail in Fig. 23 (top row), where we show an enlarged version 
of the relevant part of parameter space from Fig. 21 (right) and Fig. 22 (right). As we can 
see, the two thin extensions of the unexcluded region both lie along boundaries between areas 
where different processes have the highest statistical significance. As we have discussed, we 
would expect our method of combining channels (which only makes use of one observed limit 
for each parameter point) to be less effective at such boundaries. Thus, a dedicated analysis 
including the combination of different channels might well be able to exclude such areas. 

The results presented here update our previous result [22] in two main ways. Most 
importantly, we now use a newer version of FeynHiggs [6-10] to obtain the neutral Higgs 
self-energies, with an improved treatment of the tan /3 enhanced contributions. In addition, 
Ref. [22] uses mt = 170.9 GeV. The resulting unexcluded region B has a similar size, shape 
and position in both analyses. However, in Ref. [22], the unexcluded region A is almost 
non-existent, whereas the updated result shows a sizable unexcluded region of parameter 
space here, as we have discussed. 

Fig. 23 (top row) also shows that the bulk of the unexcluded region B lies in an area 
in which the channel — )■ bbZ (■) has the highest statistical sensitivity. The extent of 
the unexcluded region B on the higher tan /3 side is very sensitive to the h2 — )■ bb branching 
ratio, which, as we discussed in Sect. 8, is critically dependent on the /i2 — ?■ hihi decay in 
this region of CPX parameter space, since the h2 — > hihi decay width yields the dominant 
contribution to the total decay width. The extent of the unexcluded region B towards 
lower values of M/^^ is roughly determined by the edge of the region in which the channel 
h2Z —J- bbZ (■) has the highest statistical sensitivity. This boundary is also very sensitive 
to the /i2 hihi decay width, which has a large influence on the theoretical predictions of 
the other relevant channels: h2Z — t- bbZ (■), /i2/ii bbbb (■), — hihihi — t- bbbbbb ( ) 
and h2Z -> hihiZ bbbbZ (□). 

The unexcluded region B occurs within a region where the observed limit for the e^e^ — >■ 
{ha)Z — {bb)Z topology was more than one standard deviation above the expected limit 
(based on a background-only hypothesis). It is interesting to investigate the effect of this 
'slight excess' on the extent of the unexcluded region B. Fig. 24 shows what the exclusion 
would have been in the hypothetical situation in which the observed limit was exactly the 
same as the expected limit for all topologies. We see that the unexcluded region B disappears 
entirely and the size of the unexcluded region A has been reduced. We conclude that the 
presence of the 'slight excess' in the LEP results for the {ha)Z — )■ {bb)Z topology is crucial 
to the existence of substantial unexcluded regions in the CPX scenario. 

In order to further investigate the effects of our new genuine vertex corrections to the 
/i2 — hihi decay, we now compare the LEP exclusion regions based on the full result with 
cases where we have used various approximations to calculate the genuine vertex corrections 
to the /i2 — hihi decay. In Fig. 23 (centre), the genuine vertex corrections are approximated 
by the /, / corrections at zero external momenta. Fig. 23 (lower) displays the result when 
the leading Yukawa corrections to the triple Higgs vertex, as given by eqs. (76)-(89), are 
used. 

The boundary between channels related to Ifi'/j^^^P (i-e. h2Z — > hihiZ — )■ bbbbZ (□) and 
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Figure 23: The channels with the highest statistical sensitivity (left) and LEP exclusion 
regions (right) for the CPX scenario at low values of Mh^ and moderate values of tan /3. The 
top row shows a subset of the results shown in Fig. 21 and Fig. 22. In the second row, the 
genuine vertex corrections to the h2 — >■ hihi branching ratio include only the /, / corrections 
at zero external momenta and, in the third row, the vertex corrections to the /12 hihi 
branching ratio have been calculated using the Yukawa approximation. In all three rows, 
the genuine vertex corrections to the h2 — )■ hihi branching ratio are combined with the 
full propagator corrections. Left: red = hiZ — )■ bbZ [M), blue = — )■ bbZ{M), white = 
h2Z hihiZ — )■ bbbbZ (□), cyan = /i2^i bbbb (■), yellow = /i2^i hihihi — )• bbbbbb ( ), 
green = h^hi bbbb{M), purple = other channels (■). Right: green (darker grey) = LEP 
excluded, white = LEP allowed. 



h2Z bbZ (M)) and those related to Isf/i^/i^^P (i.e. /i2^i bbbb(M) and /i2^i hihihi — >■ 
bbbbbb { )) for both approximations are in about the same position as in Fig. 23 (upper). 
However, the boundaries between channels directly involving the /i2 hihi decay and those 
that do not involve this decay have shifted. In particular, the region where the channel 
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Figure 24: The LEP exclusion regions for the hypothetical case in which all of the observed 
cross section values were measured to be the same as the ones expected in a background- 
only simulation. The colour codings are: green (darker grey) = LEP excluded, white = LEP 
allowed. 

— )■ hihiZ —7- bbbbZ (□) has the highest statistical sensitivity disappears entirely if the 
'leading Yukawa' approximation is used. This has a considerable impact on the shape of the 
unexcluded region B. The /, / at = approximation, on the other hand, performs very 
well, yielding an unexcluded region B which, in position and shape, agrees very well with 
the full result. 

It is interesting to consider the effect of the experimental uncertainty on the mass of the 
top quark. Since the leading corrections to the /i2 — >■ hihi vertex are Yukawa corrections 
proportional to mf, one would expect the unexcluded region B to exhibit a strong dependence 
on rrit. The neutral Higgs masses are also very sensitive to rrit, since these Yukawa corrections 
are also the leading corrections to the Higgs self-energies. We use the experimental value 
mt = 173.1±L3 GeV [64]. Fig. 25 shows the result of running scans for rrit = 173.1 + 1.3 GeV 
and rrit = 173.1 — 1.3 GeV over the area of CPX parameter space containing the unexcluded 
region B. Parameter points which are excluded at both top masses are shown as dark green 
(darker grey), those excluded at one of the masses only are in light green (lighter grey) and 
those which remain unexcluded at both masses are in white. Therefore, the light green area 
(lighter grey area) demonstrates the uncertainty on the size and shape of the unexcluded 
region B due to the uncertainty from the measurement of the top mass. The unexcluded 
area for rrit = 173.1 + 1.3 GeV is approximately double the size of the unexcluded area for 
mt = 173.1 - 1.3 GeV. 

It is also possible to use a similar strategy to investigate the dependence of the size of 
the unexcluded region B on the choice of mt at the higher loop orders in the calculation. 
In Fig. 26, we use (a) m^"^ in the 2-loop corrections in Higgs self energies and m^"^ in the 
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Figure 25: The LEP exclusions for the CPX scenario plotted in the M/^^ — tan /3 plane. 
Results for rrit = 173.1 + 1.3 GeV and = 173.1 — 1.3 GeV. White: point is unexcluded 
in both, light green (lighter grey): point is excluded in only one, dark green (darker grey): 
point is excluded in both. 

1-loop genuine hahi,hc vertex corrections, and (b) m^^'^^ in the 2-loop corrections in Higgs 
self energies (as discussed in Sect. 5.4) and mY^'^^ in the 1-loop genuine hahbhc vertex 
corrections (as previously). Points in parameter space which are excluded for the results of 
both calculation (a) and (b) are shown in dark green (darker grey), points excluded for the 
results of one calculation only are shown in light green (lighter grey), and points in white can 
be excluded for the results of neither calculation. In this way, we can get a rough estimate of 
the dependence of unexcluded region B on higher order corrections in the t, t sector - whilst 
these corrections are significant, the existence of an unexcluded region in this part of CPX 
parameter space is confirmed in both. Note that the region in which the results of one of 
the two calculations yield an unexcluded region forks into two narrow regions at its base, 
rather than one, reflecting the different positioning of the boundaries between regions where 
different processes have the highest statistical sensitivity in the two analyses contributing to 
this exclusion plot. 

The variation of (pAt has an interesting impact on the unexcluded regions. Recall that, 
in Fig. 10 (left), we saw that there was a peak in the h2 — )■ hihi decay width at moderate 
tan /3, and a minimum at lower tan /3. Varying (p^^ by 10% has a dramatic effect on the decay 
width, through changing the magnitude and position of this peak and changing its position 
with respect to tan/3. In Fig. 27, which uses (pAt = 0.9 x n/2 and (pAt = 1-1 ^ ^/2, we can 
see these effects reflected in the /i2 — > hihi branching ratio. In particular, we see that the 
thin horizontal minimum in BR(/i2 — hihi) shifts to higher tan/3 as (pAt increases. We can 
also see a change in the shape of the region in which the /12 ^ hihi decay is kinematically 
allowed and a slight change in the lower edge of CPX parameter space as plotted in the 
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Figure 26: The LEP exclusions for the CPX scenario plotted in the M^^ —tan /3 plane. Results 
for (a) mf^ in the 2-loop corrections in Higgs self energies and mf'^ in the 1-loop genuine 
hahbhc vertex corrections, and (b) m^^'^^^ in the 2-loop corrections in Higgs self energies 

and in the 1-loop genuine hahbhc vertex corrections. White: point is unexcluded 

in both, light green (lighter grey): point is excluded in only one, dark green (darker grey): 
point is excluded in both. 



— tan /3 plane. 

As one would expect, this strong dependence of the /i2 — ?■ hihi branching ratio on (p^^ is 
reflected in a change in the balance of the processes with the highest statistical sensitivity at 
each point in parameter space as (pAt increases. In particular, the size of the region in which 
either h2Z — >■ hihiZ — >■ bbbbZ (□) or h2hi — )■ hihihi -> bbbbbb ( ) has the highest sensitivity 
decreases. The boundary between processes involving Igf^zzl"^ ^^'^ those involving Is'/j^/ij^P 
also shifts to higher tan/3. As a result, the unexcluded region B occurs at higher tan /3 as 
(pAt increases and its shape changes significantly. The unexcluded region A increases in size 
as (pAt increases. 

Fig. 28 illustrates that varying \At\ by 10% has a very substantial effect on the LEP 
exclusions in the CPX parameter space. Increasing \At\ has a large effect on the /12 hihi 
decay width. It increases the size of the peak at moderate tan /3 in Fig. 10 (left) and shifts 
the position of the minimum and the peak shown in Fig. 10 (left) to higher values of tan /3. 
Increasing \At\ also make the slope of this graph less steep above tan/3 ~ 7. We see these 
effects reflected in the h2 — )■ hihi branching ratio in Fig. 28. We also see that a higher value 
of l^tl yields a significant increase in the area of parameter space in which the /i2 hihi 
decay is kinematically allowed. 

For the case in which \At\ = 1.1 x 900 GeV, the plot illustrating the channels with 
the highest statistical sensitivity in Fig. 28 is very different from those discussed so far. 
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Colour: Br(h2->h^+h^) exclusions process with highest Stat. sens. 




Figure 27: The h2 — > hihi branching ratio, LEP exclusions and channels with the highest 
statistical sensitivity plotted in the M^^ — tan (3 plane. Upper graphs show (pAt = 0-9 x 7r/2, 
lower graphs show = 1.1 x 7r/2, other parameters taken from CPX scenario. See the 
captions of Fig. 21 and Fig. 22 for the colour codes of the plots in the second and third 
columns. 

This is partly because \g^hxh-2.z^ reduced, which drastically reduces the area where h^hx — )■ 
\M)b (■) has the highest statistical sensitivity. The area where the channel /i2^i ~^ h\h\hx — 
bbbbbb ( ) has the highest statistical sensitivity occurs at higher tan /3 than previously and is 
now unexcluded. Also, the suppression of \g'h^h.2z\'^ means that the channel h^hi — )■ bbbb{M) 
has the highest statistical sensitivity over a large region, which can only be partially excluded 
by this LEP limit. Therefore, the excluded LEP regions are dramatically different for the 
CPX scenario with \At\ = 1.1 x 900 GeV. It is worth noting, however, that this value of \At\ 
is close to an unstable region of parameter space, in which loop corrections in Higgs sector 
get extremely large. On the other hand, at \At\ = 0.9 x 900 GeV, the unexcluded region B 
has almost disappeared and the unexcluded region A has reduced in size as it is partially 
covered by two vertical excluded regions, at Mh-^ ~ 80 GeV and Mh-^ ~ 110 GeV. 

The gluino mass parameter M3 does not feature in the 1-loop corrections to the /i2 — )■ hihi 
decay or the 1-loop corrections to the Higgs masses. However, as we have seen, it can heavily 
influence the ha — )■ bb decay width. In addition, the Higgs self-energies from FeynHiggs 
depend on M3 through the 0{atas) corrections and the Arrih corrections. Therefore, it is 
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Colour: Br(h2->h^+hT) exclusions process with highest stat. sens. 




Figure 28: The h2 — > hihi branching ratio, LEP exclusions and channels with the highest 
statistical sensitivity plotted in the M/j^ — tan/3 plane. Upper graphs show \At\ = 0.9 x 
900 GeV, lower graphs show \At\ = 1.1 x 900 GeV, other parameters taken from CPX 
scenario. See the captions of Fig. 21 and Fig. 22 for the colour codes of the plots in the 
second and third columns. 

interesting to see if varying this parameter has a significant effect on the LEP exclusions. 

In Fig. 29, we show results for (pM-i = 7r/4 and (pMs = 37r/4. At = 7'"/4, the area of the 
unexcluded regions A and B increases slightly compared to the CPX setting of <^m3 = vr/2. In 
contrast, at (pMs = 37r/4, the size of the unexcluded area B decreases and the unexcluded area 
A almost disappears. Note that the parameter space is not populated above M/^^ ~ 100 GeV 
at tan/3 = 30. 

Fig. 30 illustrates the results for iMg] = 100 GeV and iMg] = 2000 GeV. Both unex- 
cluded regions are reduced in both cases. Note that, for IM3I = 100 GeV, the excluded 
region A is bounded from above. 

It is also interesting to consider the effect of varying the Higgsino mass parameter fi. The 
branching ratios shown in Fig. 31 are qualitatively very similar to those in Fig. 28. As fi 
increases, \gf^zz\'^ enhanced at the expense of |(?^\ ^P, and this determines the relative 
sizes of the regions involving these couplings. The plot with fi = 1.1 x 2000 GeV in Fig. 31 
has a large region in which the channel — hihiZ — bbbbZ (□) has the highest statistical 
sensitivity. The unexcluded regions have increased in size substantially as increases and 
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Colour: Br(h2->h^+h^) exclusions process with highest Stat. sens. 




Figure 29: The h2 — > hihi branching ratio, LEP exclusions and channels with the highest 
statistical sensitivity plotted in the M^-^ — tan /3 plane. Upper graphs show (pMs = 7r/4, lower 
graphs show (pMs = 37r/4, other parameters taken from CPX scenario. See the captions of 
Fig. 21 and Fig. 22 for the colour codes of the plots in the second and third columns. 

they have joined up. 

Similarly, the effect of decreasing the soft-breaking term Msusy by 10%, as shown in 
Fig. 32, can be explained by an enhancement of Igf^zzl"^ expense of \g'^^hiz\'^ ^ 

suppression of BR(/i2 — )■ bb) as Msusy decreases. This results in a very large unexcluded 
region, covering the majority of the parameter space above tan /3 ~ 9. In contrast, increasing 
the value of Msusy to Msusy = 1.1 x 500 GeV greatly reduces the amount of unexcluded 
parameter space. In addition, note that the parameter space is populated at lower values of 
tan/3 at this value of Msusy- 

10.3 Comparison with CPsuperH 

As we discussed in Sect. 5, when comparing our calculation with the results from CPsuperH, 
we need to make sure that we account for the difference between the definitions of the input 
parameters. 

When comparing LEP exclusion plots, we shall investigate both the CPX and CPX 
scenarios, as defined in Sect. 3 (recall that, in the CPX scenario, the input parameters At, 
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Colour: Br(h2->h^+hT) exclusions process with highest stat. sens. 




Colour: Br(h2-^h^+h^) exclusions process with highest stat. sens. 




Figure 30: The h2 — > hihi branching ratio, LEP exclusions and channels with the highest 
statistical sensitivity plotted in the M/j^ — tan /3 plane. Upper graphs show IM3I = 100 GeV, 
lower graphs show IM3I = 2000 GeV, other parameters taken from CPX scenario. See the 
captions of Fig. 21 and Fig. 22 for the colour codes of the plots in the second and third 
columns. 

Ml and are defined in the DR scheme, while in the CPX scenario, they are defined 
in the on-shell scheme). The CPX scenario is used most frequently in the literature as 
an example of an interesting CP-violating MSSM scenario and closely resembles the CP- 
violating scenario used by the LEP Higgs Working Group [2]. The CPX scenario is a more 
natural scenario in which to perform a Feynman-diagrammatic calculation and therefore has 
been used throughout the majority of this paper, up to this section. 

However, we shall begin the discussion by an investigation of two scenarios with less 
extreme parameter values than the CPX and CPX scenarios, i.e. where the higher-order 
effects are expected to be smaller. We will use the (CP-conserving) m'^^^ scenario, where 
we interpret the stop sector parameters At, Ml and M^j^ as DR parameters (tan/3 = 10, 
Mh± = Mp^(Ms) = M^{Ms) = 1000 GeV, fi^{Ms) = 200 GeV, M3 = 800 GeV). We 
also consider an 'm™^^-like' scenario with CP violation (as for m™'^'^, but ^^f" — 0a/^ — ^/2)- 
Since we have defined At, M|, M?^ and /i in the DR renormalisation scheme at Ms, we must 
make use of the parameter conversions described in Sect. 5, in order to obtain yi°ii-sheii^ 
j^^on-sheii^ j^^on-sheii ^^^^ /i°^(m(), which are the input parameters that need to be given to 
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Colour: Br(h2->h^+h^) exclusions process with highest stat. sens. 




Figure 31: The h2 — > hihi branching ratio, LEP exclusions and channels with the highest 
statistical sensitivity plotted in the M/^^ — tan/3 plane. Upper graphs show fi = 0.9 x 
2000 GeV, lower graphs show fi = 1.1 x 2000 GeV, other parameters taken from CPX 
scenario. See the captions of Fig. 21 and Fig. 22 for the colour codes of the plots in the 
second and third columns. 

FeynHiggs. 

In Fig. 33 (a), we display the result from CPsuperH (blue dashed) and the result from 
FeynHiggs using 0{as) and 0{at) contributions to the parameter shifts (black solid) (as 
discussed in Sect. 5) in the m™^^ scenario. We can see that the maximum value of Mh-^ 
agrees to within 0.1 GeV between the two codes, although the maximum occurs at a value 
of Af^ which is ~ 200 GeV higher for the FeynHiggs result than the CPsuperH result. Note 
also that there is a somewhat larger discrepancy between the two codes at small values of 
Af^, where CPsuperH predicts a mass ~ 3 GeV higher than FeynHiggs. Since FeynHiggs 
contains the full 0{at<ys), 0{af) and the complete 1-loop corrections to the neutral Higgs 
self-energies in the MSSM with real parameters, while CPsuperH only contains the leading 
logarithmic contributions, we should not expect our 0{as) and 0{at) parameter shifts to 
convert perfectly between the two calculations. 

Also shown in Fig. 33 (a) is the result where the FeynHiggs calculation is parametrised 
in terms of the MS top mass (magenta dotted), as discussed in Sect. 5.4, instead of the 
on-shell top mass (black solid) . This changes the Higgs mass by less than 2 GeV (for Af^ < 
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Figure 32: The h2 — > hihi branching ratio, LEP exclusions and channels with the highest 
statistical sensitivity plotted in the M^^ — tan/3 plane. Upper graphs show MgusY = 0.9 x 
500 GeV, lower graphs show MgusY = 1.1 x 500 GeV, other parameters taken from CPX 
scenario. See the captions of Fig. 21 and Fig. 22 for the colour codes of the plots in the 
second and third columns. 

3600 GeV), which is appropriately small since this is formally a higher order effect. 

Another contribution to the remaining discrepancy between the results for the Higgs 
mass predictions in the two codes is the difference in the way that the electric charge is 
parameterised at lowest order in the two codes. We investigate the magnitude of this effect 
in the following way. The input to CPsuperH is supposed to be a{Mz) (blue dashed), but we 
can investigate the effect of changing it to a —t- GF'2.^/2M^rs'^r/{4:^^) (green dotted), which 
should only affect the calculation at the level of the unknown higher order contributions. In 
Fig. 33 (b), we can see that this changes the lightest Higgs mass by less than 1.3 GeV for 
Af^ < 3600 GeV. 

A meaningful comparison between the two Higgs mass calculations can also be carried 
out in the m™*^^ scenario by neglecting the parameter shifts in Mi and M^^ and using the 
values of X^^ and X^^, which maximise the lightest Higgs masses from FeynHiggs and 
CPsuperH, respectively (corresponding to independently defining a m™*^^ scenario in the on- 
shell scheme and a m™*^^ scenario in the DR scheme). Using the leading 1-loop and 2- loop 
corrections, it is possible to predict that this will occur at approximately X^'^ = 2000 GeV 
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(a) (b) (c) 

Figure 33: Lightest neutral Higgs mass as a function of Af^ in the m^^^ scenario, with 
tan/3 = 10, Mh± = 1000 GeV. All FeynHiggs results shown here use the full Higgs self- 
energies, and the DR to on-shell parameter conversion involves terms at both 0{as) and 
and the calculation is parametrised in terms of m^^ unless otherwise stated. 

(a) blue dashed: CPsuperH, black solid: FeynHiggs, magenta dotted: FeynHiggs with code 
edited such that calculation is parametrised in terms of rut. 

(b) blue dashed: CPsuperH, green dotted: CPsuperH but with the result parametrised in 
terms of the Fermi constant instead of a{Mz), black solid: FeynHiggs. 

(c) blue dashed: CPsuperH as a function of Af^ as before, red dotted: FeynHiggs as a 
function of Af^. 

and X^=y6x 1000 GeV, see Ref. [77]. In Fig. 33 (c), we plot the CPsuperH result for Mh, 
against Af^ (blue dashed, as before), together with the FeynHiggs result for M/,,^ against Af^ 
(red dotted). The maxima in Mh, occur at the points (Xf ^, MhJ = (1990 GeV, 129.3 GeV) 
and (Xp,M/,J = (2350 GeV, 129.2 GeV), i.e. this definition of the m'j^'''' scenario does 
indeed give remarkable agreement between the two codes for the value of Mh-^. The two 
vertical lines show Xj^^ = 2000 GeV and = \/Q x 1000 GeV (we can see that using 
these approximate values of X^^ and X^^ predicts lightest Higgs masses which are within 
0.1 GeV of the true values at the maxima). 

It is interesting to compare the corresponding shift in At i.e. AAt = AXt = X^^—X^^ = 
2350 — 1990 = 360 GeV, with the shift AAt which we obtain from performing the parameter 
conversion. At X^^ = \/Q x 1000 GeV, AAt = 510 GeV when just 0{as) corrections are 
included in the parameter conversion and AAt = 600 GeV when both 0{as) and 0{at) 
corrections are included. Comparison with the analysis in Ref. [77], which was performed 
for the leading 1-loop and 2-loop contributions, shows a similar behaviour if only the 0{as) 
correction to the parameter conversion is taken into account. Incorporating also the 0{at) 
correction turns out to give rise to a slightly larger discrepancy between the shifts AAt and 
AAt- The numerical effect on the lightest Higgs mass, however, is small over this range of 
At. 

As discussed in Sect. 4, the only 2-loop pieces contained within FeynHiggs which have 
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Figure 34: Lightest neutral Higgs mass as a function of Af^ in the m™^^ scenario (left) and a 
CP-violating m™^''-like scenario (right), where (p^-m = (pMs = 7r/2. In both plots, tan/3 = 10 
and Mh± = 1000 GeV. All FeynHiggs results use Higgs self-energies up to 0{atas) only 
(with the exception of the orange (thin solid) line) and the calculation is parametrised in 
terms of mf^. Blue dashed: CP super H as in Fig. 33, red dotted: FeynHiggs result without 
parameter shifts for the input parameters i.e. Ap = 0, green solid (thick): FeynHiggs with Ap 
at 0{as), purple dot-dashed: FeynHiggs result with Ap involving contributions at 0{as) and 
0{at), orange solid (thin): as for purple, except subleading 2-loop corrections are included 
in the Higgs self-energies, calculated at all possible combinations of (pxt = 0, vr, (px^ = 0, tt 
and (f>M3 = 0, TT. Therefore, there are eight orange lines in total, but they lie too close to 
each other to be distinguishable by eye. 

been calculated with explicit phase dependence are of 0{atas). Before we start to discuss 
the CP-violating MSSM, it is instructive to compare the results from CPsuperH with the 
results given by FeynHiggs when the only 2-loop pieces which are included are 0{atas)- 
This is shown in Fig. 34 (left). Also shown are the FeynHiggs results when the parameter 
shifts are neglected (red dotted), in order to give an idea of the size of the parameter shifts. 
We have a choice for the parameter shifts: we can either use the shifts at both 0{as) 
and 0{at) (purple dot-dashed) as before, or restrict to 0{as) pieces (green solid (thick)). 
Theoretically, restricting to 0{as) pieces (green solid (thick)) is preferable, since neither 
Higgs mass calculation now involves the full set of 0{a^) corrections. It can be seen that 
this choice makes very little difference in this scenario. We can see that, when we restrict 
to 2-loop pieces at 0{atas) in FeynHiggs, the resulting Higgs mass at the peak in the plot 
is approximately 7 GeV lighter than when all available corrections in FeynHiggs are used 
(and therefore also significantly lighter than the mass given by CPsuperH oX the peak). In 
Fig. 34 (right), we introduce CP violation: this figure has the same parameters as Fig. 34 
(left) (tan/3 = 10, Mh± = M^{Ms) = Mff{Ms) = 1000 GeV, fi^iMs) = 200 GeV, 
IM3I = 800 GeV), except that At, A^ and M3 are now complex: 0^dr = (f)At = = 7r/2. 
Note that both the 0{as) and 0{at) parameter shifts include the full phase dependence, as 
described in Sect. 5. We can see that, in this scenario, the introduction of complex phases 
does not significantly affect the lightest Higgs mass. Therefore, we can deduce that the 
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dominant source of the discrepancy between the masses in this CP-violating example are the 
neglected 2-loop corrections. It is instructive to consider whether the calculation of these 2- 
loop corrections to the neutral Higgs self-energies for the MSSM with real parameters can be 
used to give an estimate of these corrections in the complex MSSM. It is possible to evaluate 
these additional corrections at a phase of and vr for each complex parameter and then 
perform an interpolation between these results (this is an optional feature in FeynHiggs). 
In Fig. 34 (right), we show predictions for the lightest Higgs mass (orange solid (thin)) in 
the CPX scenario when the neutral Higgs self-energies are calculated up to 0{atas) with 
full phase dependence (as previously) and the subleading 2-loop corrections are calculated 
at combinations of (pxt = 0, vr, (px,, = 0, tt and (pMs = 0, n. These eight combinations give 
very similar results, and their curves can not be distinguished in the plot, and, obviously, 
an interpolation between these results would also lie on this same line. It is, however, 
difficult to evaluate the accuracy of this Higgs mass prediction. On the one hand, it provides 
better agreement with the CPsuperH lesult, and the fact that the eight CP-conserving limits 
give very similar results could be interpreted as an indication that a fully phase-dependent 
calculation of these corrections may have only a mild dependence on the phase. However, a 
confirmation of this interpretation would require the actual result of a fully phase-dependent 
calculation. In an 'extreme' scenario like the CPX scenario, which is close to unstable regions 
of parameter space, the evaluation of the combinations of real parameters needed as input 
for the interpolation can be problematic. In particular, higher-order contributions in the 
CP-conserving limit (pxt = turn out to be unphysically large in the CPX scenario and 
cause the mass calculation to fail (to obtain a real and positive mass value, the limit (p^^ = 
must be used instead). Since we will now turn our attention to such 'extreme' examples, 
exhibiting very large higher order corrections, we shall restrict in these cases the Higgs self- 
energies calculation in the complex MSSM to 0{atas), as we have done throughout the rest 
of the paper. 

We shall first consider a CP-conserving version of the CPX scenario i.e. the CPX sce- 
nario with all phases set to zero (and fi^^{Ms) = 2000 GeV), as this will exhibit the main 
characteristics of the CPX scenario (and, since we have more available corrections in the 
MSSM with real parameters than the MSSM with complex parameters, it will be easier for 
us to investigate). In Fig. 35, we compare the CPsuperH result (blue dashed) to the full 
FeynHiggs result with Ap at 0{as) and 0{at) (black solid (thin)), the FeynHiggs result with 
Higgs self-energies only up to 0{atas) and Ap at 0{as) and 0{at) (purple dot-dashed), the 
FeynHiggs result with Higgs self-energies up to 0{atas) and Ap only at 0{as) (green solid 
(thick)). Note that, when the 0{at) contributions to Ap are included, this incorporates 
the conversion of fjP^{Ms) to fj,^^{mt), using eq. (58). Also illustrated is the FeynHiggs 
result when the parameter shifts are neglected (red dotted), in order to demonstrate the 
numerical effect of the parameter shifts. In each plot in Fig. 35, we can see that the effect 
of the parameter shifts at 0{at) on the mass calculation becomes very large once the nom- 
inal values of this CP-conserving CPX-like scenario are approached (Ml'^.jy = 500 GeV, 
= 2000 GeV, Af^ = 1000 GeV, Mh± = 140 GeV). Indeed, for the nominal values of 
this scenario, the lightest on-shell stop mass is driven to imaginary values when both the 
0{as) and 0{at) parameter shifts are used (purple dot-dashed and black solid (thin)). These 
unphysically large effects are an indication that this scenario is close to an unstable region, 
which causes instabilities in the parameter conversion (this could imply that it would be 
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Figure 35: Lightest neutral Higgs mass in a CP-conserving version of the CPX scenario at tan /3 = 10, Mu± = 140 GeV, as 
function of (left to right) M^^^gyi^s), fJ'^^iMs), \Af^{Ms)\ and Mh±. Blue dashed: CPsuperH result, red dotted: FeynHiggs 
with Higgs self-energies up to 0{atas), with Ap = 0, green solid (thick): FeynHiggs with Higgs self-energies up to 0{atas) 
g corrections, with Ap at 0{as), purple dot-dashed: FeynHiggs with Higgs self-energies up to 0{atas) corrections, with Ap at 

0{as) and 0{at), black sohd (thin): FeynHiggs with full Higgs self-energies, with Ap at 0{as) and 0{at). 
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Figure 36: Lightest neutral Higgs mass in the CPX scenario at tan /3 = 10, Mh± = 140 GeV, as function of (left to right) 
M^^gyi^s)-, A^; l^?^(^5)l ^iid Mh±. (See Fig. 35 for the definitions of the abbreviations labelling the different curves.) 



appropriate to explicitly decouple the heavier particles in the conversion to the DR scheme). 
In order to avoid these instabilities, we will, in the following, restrict the parameter conver- 
sion to 0{as)- Application of this parameter conversion, in conjunction with 2-loop mass 
corrections at 0{atas) (green solid (thick)), yields a relatively good agreement between the 
FeynHiggs and CPsuperH results. 

We observe the same quantitative behaviour when we switch on the CP-violating phases 
in this scenario (Fig. 36), and therefore we can use the knowledge gained from the CP- 
conserving version to help interpret these characteristics. Once again, we find that the 
parameter conversions at 0{as) and 0{at) result in unphysically large corrections to the 
lightest Higgs mass and therefore restrict the parameter conversion to 0{as) (green solid), 
which compares relatively well to the CPsuperH result (blue dashed). 

We will now carry out a complementary analysis in the CPX scenario, which is defined 
in terms of on-shell stop sector parameters. Therefore, the parameter shifts are now used 
to convert from on-shell to DR parameters, which can then be given as input to CPsuperH. 
Results for the lightest Higgs mass are shown in Fig. 37 for CPsuperH with no parameter 
conversion (blue dashed) and CPsuperH with parameter conversion at 0{as) (green dot- 
dashed). Also shown is the FeynHiggs result where the masses are calculated at C(ajas) (red 
dotted). Once again, we can see that this scenario lies near the border of stable parameter 
space, which is necessary if we would like to investigate very low lightest Higgs masses. We 
find that, when the parameter conversion is performed in the on-shell to DR direction, the 
effect of the parameter shifts tends to push the results into more stable regions of parameter 
space (rather than into more unstable regions, as it happens in the case of the CPX scenario). 
However, we have confirmed that this parameter conversion still gives rise to large stop mass 
corrections at 0{at). 

Before starting a comparison of the LEP exclusions in CPX or CPX parameter space, it is 
useful to describe the method that was used by the LEP Higgs Working Group to tackle this 
issue of consistently comparing results calculated in using different renormalisation schemes. 
The LHWG analysis was carried out in the CPX scenario and both FeynHiggs version 2.0 [10] 
and CPH [5] (a predecessor of the program CPsuperH [11,12]) were considered. The relation 
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was used [132] to convert between different definitions for \At\, with = 0.108 and Xt = 
\At\ — nl tau/S. ^ The shifts in (j)At and Msusy(= Ml = Mi^) were neglected, and Ai, was 
set to be the same as At. 

Near the region Mh^ ~ 45 GeV, tan /3 ~ 7, this approximation gives a larger shift in 
\At\ than would be given from the full 0(«s) shifts, and therefore increases the agreement 
between the two Higgs mass codes. The LHWG then carried out two parallel analyses, 
one using Higgs sector predictions from CPH and one using Higgs sector predictions from 
FeynHiggs (apart from the triple Higgs coupling, as discussed in Sect. 10). Since they found 
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Ms Mf Ms V^/. 



MS 
t 



given by [77]: X 



was obtained for nig 



8- 



, where Ms = 



SUSY "H I 
2 



Ms and the assumptions mt/Ms ^ 1 and mtXt/Mg <C 1. 
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Figure 37: Lightest neutral Higgs mass in the CPX scenario at tan/3 = 10, Mh± = 140 GeV, as function of (left to right) M^^gy, 
fj,, \Af^\ and Mh±. Red dotted: FeynHiggs with Higgs self-energies up to 0{atas) corrections, blue dashed: CP super H with 
Ap = 0, green dot-dashed: CP super H with Ap at 0{as), 
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Figure 38: The /12 h\hi branching ratio, LEP exclusions and channels with the highest 
statistical sensitivity for the CPX scenario, plotted in the M^^ — tan /3 plane. Top row: 
Using the calculations presented in this paper, parameter conversions at 0{as) as described 
by eqs. (48)-(53), and neutral Higgs self-energies from FeynHiggs up to 0{atas)- Bottom 
row: Using Higgs masses, couplings and branching ratios calculated by CPsuperH [11, 12]. 
See the captions of Fig. 21 and Fig. 22 for the colour codes of the plots in the second and 
third columns. 



significant differences between the two results, they considered a point in CPX parameter 
space to be excluded only if it could be excluded in both the analysis based on the CPH 
results and the analysis based on the FeynHiggs results. 

In Fig. 38 (top row), we use the topological cross section limits from LEP in conjunction 
with updated results for the Higgs masses, couplings and branching ratios, in the CPX 
scenario using the parameter conversions at 0{as). In particular, we use our complete 1- 
loop diagrammatic calculation for the hi — )■ hjhk decay processes with full phase dependence 
as described in Sect. 6 and we use renormalised neutral Higgs self-energies obtained from 
FeynHiggs (which includes corrections at 0{atas) with full phase dependence). We compare 
this to the results using the Higgs masses, couplings and branching ratios from the program 
CPsuperH, shown in Fig. 38 (bottom row). The CPsuperH prediction for the h2 — )■ hihi 
branching ratio is qualitatively similar to our result, although the thin 'knife-edge' region is 
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Colour: Br(h2^h^+h-|) exclusions process with highest stat. sens. 




Figure 39: The h2 — ?■ hihi branching ratio, LEP exclusions and channels with the highest 
statistical sensitivity for the CPX scenario, plotted in the M/^^ — tan /3 plane. Top row: Using 
the calculations presented in this paper and neutral Higgs self-energies from FeynHiggs up 
to 0{atas)- Bottom row: Using Higgs masses, couplings and branching ratios calculated by 
CPsuperH [11,12] and parameter conversions at 0{as) as described by eqs. (48)-(53). See 
the captions of Fig. 21 and Fig. 22 for the colour codes of the plots in the second and third 
columns. 

at tan /3 ~ 3 rather than tan /3 ~ 4.5 and the peak at moderate tan /3 is higher but narrower, 
in particular, around M/j^ ~ 45 GeV. Also, different processes have the highest statistical 
sensitivity at different parts of the parameter space (Fig. 38 right column). Most significantly, 
the area where the process — bbZ (■) has the highest statistical sensitivity is larger 
when using the Higgs sector predictions outlined in this paper, rather than the predictions 
given by CPsuperH. As a result, the unexcluded part of parameter space is much larger 
for M/jj < 60 GeV when using our calculation and the unexcluded regions join together. 
However, for Mh^ > 60 GeV, the unexcluded area is larger when the CPsuperH results are 
used. 

Finally, in Fig. 39, we compare the results from CPsuperH and our calculation (as shown 
previously in Fig. 16 (left). Fig. 21 (right) and Fig. 22 (right)) for the Higgs cascade decay 
width and LEP exclusions in the CPX scenario (defined in terms of on-shell stop sector pa- 
rameters). Again, the CPsuperH piediction for the /12 hihi branching ratio is qualitatively 
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similar to our result, although, once again, the thin 'knife-edge' region is at lower tan /3 and 
the maximum at moderate tan/3 is higher. As we saw in Fig. 38, the CPsuperH result has a 
large region where hiZ — )■ bbZ (■) has the highest statistical sensitivity at medium-to-high 
tan /3 and M/^^, which is not seen in our result. Once again, the area where processes involv- 
ing /11/12 pair-production (■ and ) have the highest statistical sensitivity is larger in the 
CPsuperH result than in ours. As a result, the region where — > bbZ (■) has the highest 
statistical sensitivity does not significantly overlap with the region with large /12 ~^ hihi 
branching ratio, and therefore the region M/,^ < 65 GeV is almost entirely excluded. How- 
ever, once again, the CPsuperH result has a larger unexcluded region at medium-to-high 
tan /3 and Mh^ . 

In summary, we have compared the results presented in this paper with those given 
by the program CPsuperH. Since these results are very sensitive to the Higgs masses, we 
started with a comparison between the masses given by FeynHiggs and CPsuperH (including 
appropriate parameter transformations, which we have extended to the complex MSSM) 
in the m^^^ scenario, in a m™^^-like scenario with added CP violation, the CPX scenario 
(in which the stop sector parameters have DR definitions), the CPX scenario with vanishing 
phases, and the CPX scenario (in which the stop sector parameters have on-shell definitions). 
We find good agreement between the mass calculations in the rn^^^ and m^'^^-like scenarios 
once the appropriate parameter shifts are taken into account. The CPX scenario has been 
used by the LEP Higgs Working Group in their analysis and is the most frequently used 
CP- violating scenario in the literature. We find that the 0{at) corrections in the parameter 
transformations for the stop masses are very large in this scenario and drive the ti mass to 
unphysical values. We interpret this as an indication that this scenario is close to an unstable 
region of parameter space (i.e. a region with relatively large higher-order corrections), where 
we must be careful to avoid unphysically large corrections in the parameter shifts (which 
may require to explicitly decouple heavy particles in the conversion). When restricting to the 
0{as) parameter conversions, we find rough agreement between the mass predictions, and we 
observed similar behaviour when we considered the CPX scenario with appropriate parameter 
conversions at 0{as)- The variation in the Higgs sector predictions given by our calculation 
(using neutral Higgs self-energies from FeynHiggs) and CPsuperH are refiected in significant 
differences in the size and shape of the regions of parameter space that could not be excluded 
by the LEP results. When the topological cross section limits from LEP are confronted with 
the predictions of our calculation, significantly less parameter space can be excluded at low 
Mh^. However, there is qualitative agreement about the existence of an unexcluded region 
at ~ 45 GeV. Using either our calculation or the CPsuperH calculation results in 
large unexcluded regions at medium-to-high tan /3 and M/^ . In general, more of the CPX 
parameter space can be excluded by LEP, as compared to the CPX parameter space. Once 
again, both calculations are consistent with an unexcluded region at M^^ ~ 45 GeV in the 
CPX scenario, although, when CPsuperH Higgs sector predictions are used, this region is very 
small, and might be covered by a dedicated analysis of this scenario which could make use of 
a combination of different LEP Higgs search channels (as discussed in Sect. 10.2). However, 
both calculations once again yield large unexcluded regions of parameter space at medium- 
to-high tan /3 and M/^^ . Since the shape, size and position of the unexcluded regions vary 
significantly depending on which method was used to calculate the Higgs sector predictions, 
we adopt the LEP Higgs Working Group philosophy and only consider a point in parameter 



65 



CPX 



CPX 




Figure 40: Combining results for the LEP exclusions in the CPX (left, using Fig. 38) and 
CPX (right, using Fig. 39) parameter space. Green (darker grey): point in parameter space 
can be excluded using the topological LEP limits and the Higgs sector predictions from both 
our calculation (which uses neutral Higgs self-energies from FeynHiggs) and the calculation 
implemented into CPsuperH, white: is excluded by only one or by neither of the calculations. 

space excluded if it can be excluded using both methods independently. We therefore show 
these 'combined' exclusion plots in Fig. 40 for the CPX and CPX scenarios. We conclude 
that there is still a sizeable area of parameter space in these scenarios below the kinematical 
limit at LEP which the LEP Higgs topological cross section limits are unable to exclude (in 
agreement with the LEP Higgs Working Group analysis) and therefore these scenarios with 
a very light neutral MSSM Higgs boson remain phenomenologically very interesting. 



11 Conclusions 

In this paper we have presented theoretical predictions concerning the Higgs sector in the 
MSSM with complex parameters, which are especially relevant when considering MSSM 
scenarios that are particularly challenging for the past and present Higgs searches at LEP, 
the Tevatron and the LHC. 

We have calculated the full 1-loop vertex corrections within the Feynman-diagrammatic 
approach for the process ha — )■ h^hc, taking into account the dependence on all complex 
phases of the supersymmetric parameters and the full momentum dependence. We have 
included the full propagator corrections, using neutral Higgs self-energies as provided by 
FeynHiggs [6-10], and we have consistently included 1-loop mixing with the Z boson and 
the unphysical Goldstone-boson degree of freedom. Our results are currently the most precise 
predictions for the ha — )■ hbhc decay width. These results will be included in the publicly 
available program FeynHiggs. 
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We have found that the genuine vertex corrections to the triple Higgs vertex are numer- 
ically very important. This holds both for the cases of CP conservation and CP violation. 
The inclusion of the genuine vertex corrections changes the predictions for the decay widths 
drastically as compared to an approximation based solely on propagator-type corrections. 
Using the leading Yukawa contributions or the full fermion/sfermion vertex corrections at 
zero incoming momentum yields a prediction for the decay width which is closer to the full 
result, but we still find deviations of ~ 30% and ~ 15%, respectively, in the example of 
the CPX parameter space at tan/3 ~ 9, Mh^ ~ 40 GeV. Yet, these approximations can 
be useful in terms of effective couplings. The Yukawa approximation including all leading 
corrections can be expressed in a very compact form, thus providing a convenient way to 
go beyond the tree level vertex contributions. The effective coupling created from the full 
fermion/sfermion contributions at vanishing external momentum is a more sophisticated 
effective coupling approximation. These effective couplings can be used for theoretical pre- 
dictions of cross sections for processes such as e"^e~ — hiZ — hihiZ at a future Linear 
Collider and the corresponding process at the LHC, which are of great interest in order to 
directly access the Higgs self-couplings and thus investigate a crucial element of the Higgs 
mechanism. We find that the genuine vertex corrections enhance the triple-Higgs squared 
coupling of the light CP-even Higgs boson of the MSSM by 20-60% in the phenomenologi- 
cally relevant region as compared to the lowest order coupling in the Standard Model for the 
same value of the Higgs mass. This could have interesting implications for the sensitivity of 
searches at the LHC and a future Linear Collider to effects of the triple-Higgs coupling. 

We also have presented the full 1-loop electroweak vertex corrections to the ha ff de- 
cay width in the complex MSSM, including full phase dependence. We have supplemented 
these contributions with 1-loop QED and SM QCD corrections, resummed SUSY QCD 
contributions, propagator corrections calculated using neutral Higgs self-energies from Feyn- 
Higgs, and 1-loop propagator mixing with Goldstone bosons and Z bosons. Our calculation 
has already been incorporated into the program FeynHiggs. 

We performed on-shell to DR parameter conversions at 0{as) and 0{at) for the general 
case of complex supersymmetric parameters, to be used when comparing results of calcula- 
tions based on different renormalisation schemes. We also present a simple approximation 
to the 0{as) contribution to these parameter shifts, which is valid in the complex MSSM. 

We investigated the numerical impact of parametrising the neutral Higgs self-energies (in 
the Feynman-diagrammatic approach) in terms of the MS top mass, rather than the on-shell 
top mass, which is formally a 3-loop effect. This will be incorporated as a possible option in 
FeynHiggs. 

Furthermore, the investigated t, t, b, b genuine vertex corrections to the LEP Higgsstrahlung 
and LEP pair production processes in the CP-violating MSSM were found to have a small 
numerical impact in the CPX parameter space. 

Using our new theoretical predictions in conjunction with the topological cross section 
limits from the LEP Higgs searches (there exist no relevant limits for this parameter re- 
gion from the Tevatron), we were able to investigate the effect of the new contributions on 
the exclusion regions in the CPX benchmark scenario. Of particular interest is the region 
30 GeV ^ M/jj < 50 GeV, 3 ^ tan/3 ^ 10, which previous analyses had not been able to 
exclude, despite the relatively low values of Mh^ involved. Since the coupling of the lightest 
Higgs boson to gauge bosons is suppressed in much of this region, processes involving the 
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second heaviest or the heaviest neutral Higgs boson are important. Across the majority of 
this region, the /i2 — hihi decay is dominant and therefore a precise theoretical prediction 
for this decay width is crucial for investigating whether an unexcluded region actually exists 
and for mapping out its extent. Our results confirm the existence of an unexcluded region at 
M/ij ~ 45 GeV. Furthermore, we find that there is a rather large unexcluded region in the 
CPX-type scenarios below the LEP limit for a SM-like Higgs, occurring at medium to high 
values of tan /3 and extending down to M/^^ values of about 60 GeV. We have performed a 
detailed analysis of the dependence of the unexcluded parameter regions on variations of the 
MSSM parameters. 

After discussing in detail the issues involved in such a comparison, we have compared 
our results for the impact of the existing Higgs search results on the parameter region of 
the CPX scenario with the results obtained using the publicly available program CPsuperH. 
While both analyses confirm the existence of unexcluded regions in similar parts of the CPX 
parameter space, sizable differences occur in the detailed extent of these regions. 

Our results should serve as a motivation to focus Higgs searches at the LHC not only 
on the mass range above about 114 GeV, corresponding to the allowed region for a SM-like 
Higgs, but to perform dedicated searches also for much lighter Higgs states. 
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A Relations between T^h^G ^ind H/^.^ 

Using a Feynman-diagrammatic calculation, we find that the following relations hold at 
1-loop order in the -R^-gauge: 

MzT.^,g{P^) + Zp^S^j^(p^) = — ^ S,3tA, 

MzT.AG{p') + ip'^Az{p') + Mz{p' -m\)U{p') = —^—{cpt^^-spt^,), 

2 

Mw^H-G+{.P^) - P^^H-w+{.P^) + Mw{p^ - m]j±)f^{p^) = -^{ci^t^^- spt^^+itA) , 

2 

MwT.G-G^{p')-'^p'^G~W+{p')-Yl-^Ww{p') = 77^M</.2+C/3t^i), (127) 
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where 

foip') = - .^ ? ,.2 Sfi-aC^-aMUz [^0 {p\ ^1 M^z) " B, {p\ mj,, Mliz)] , (128) 



/o, /± are finite and do not contribute if tlie particle is on-sliell. 
In terms of renormalised quantities, tlie relations are given by 

2^ , / 2x 



(129) 



^ag{p') + ^^Azip') + ip' - m\)Up') 



^gg{p') + '^-^taz{p') - ^tUp') 
z ^ 



^H-G+{P^) - T^^H-W+{P') + {P' - m\j±)f±{j)' 



^G~G^{p') - ^f^G-W^{p') - ^SW(P') 



0, 


(130) 


0, 


(131) 


0, 


(132) 


0, 


(133) 


0, 


(134) 


0. 


(135) 



These relations have been checked algebraically using the t, t, 6, h sector and the gauge and 
Higgs boson sector. These relations have also been checked numerically for the entire MSSM. 
These expressions first appeared in Ref. [133] and were derived using a BRST transformation 
in Ref. [134]. 

B Explicit form of counterterms 

In this section, we briefly summarise the counterterms used in this paper. 
We renormalise the gauge bosons on-shell: 



5Ml = ReSL(M|), (136) 
6M^ = ReS^^(M^), (137) 
6Zvv = -ReS'^^, (138) 

SZ,z = -2 (139) 

SZz, = 2^^, (140) 

where 'T' denotes the transverse part of the self-energy, and as before the prime indicates 
the derivative w.r.t. the external momentum squared. We use a combination of on-shell and 
DR renormalisation for the Higgs sector, as in Ref. [7]: 
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5M^± = ReSH±(M^±), (141) 
6Zu, = -[ReS;^^j'^\ (142) 
SZn. = -[Re^^X^. (143) 

'div' indicates that we have just kept the terms proportional to the divergent piece — 
'Je + log(47r). We define the tan/3 counterterm using 

tan/3 tan/3(l + 5tan/3) =tan/3(l + ^(5Z^2 -5Z^J). (144) 

This choice has been shown in Refs. [72, 113, 135] to yield numerically stable results. We also 
need to fix the renormalisation scale for 6Z-}i^,6Zy^^, which we choose to be fij-^n = ^t, as in 
Ref. [7]. We set 

STh. = -T,„ (145) 

to ensure that the renormalised tadpoles vanish. 

In all of our calculations involving the full complex MSSM, we choose to parametrise the 
result in terms of a{M^), where a(M|) = a{0)/ (1 — Aa), and Aa is a finite quantity that 
can be split into the contribution from the e, /i, r leptons and the contribution from the light 
quarks (i.e. all quarks except t), Aa = Aaicpt + Aa[^j. Aaicpt has been calculated to 3-loop 
order [136] as 

Aaiept = 0.031497687, (146) 
while Aa^^^ has been determined using experimental data via a dispersion relation [137] as 

Aa^^^^ = 0.02755 ± 0.0023. (147) 
We therefore use the charge counterterm: 



7e(M|) _ v^allloops^„2N| ^ v light /in loops/ ^2 \ 

~ 2 9^2^77 y1 )\q2=0 2 ^77 



+ 2i^^"^"''''"'°°"(^^)+S^Mr' ^^^^^ 



Note that the second term in the right-hand side of eq. (148) cancels the light fermion con- 



"light/ 



rcn 



tributions in the first term, which involve large contributions proportional to a log 

arising from the running of a from = to a higher energy scale, which are problematic 
because the masses of the light quarks are not well defined. 

For calculations which only involve Standard Model fermions and their superpartners in 
loops, we choose to parametrise the electric charge in terms of the Fermi constant Gp (as 
used, for example, in Ref. [138]), yielding the charge renormalisation counterterm 

SZ^^ = ^-T^ - ^ i^lwiO) - SM^w) ■ (149) 

SW ^My^ 
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We renormalise the quark sector on-shell using the conditions^ 



6m 
6Z^ 



-Re [m (S^(m2) + S^(m')) + SV') + S''(m')] 



Re 



1 
2m 



m 



5Z 



R 



Re 



2m 



-m" ( S^' (m") + S^' (m") ) - m ( S'' (m") + S"' (m' 



(150) 



(151) 



(152) 



In order to perform the DR to on-shell parameter conversion in Sect. 5, we apply the 
renormalisation conditions whi 
rections to the Higgs self-enerj 
mggs [6-10]: 



22 



12 



21 



were 
[110] 


used in the calculation of the leading 
, which were incorporated into the pro 


2-loop cor- 
gram Feyn- 


ti 


= ReEi^^(m|), 


(153) 




= ReEi^2(m|), 


(154) 




= l(ReS,-^,(m|) + R;S,^,(m|)), 


(155) 






(156) 



C Triple Higgs vertex corrections and Higgs self-energies 
in the Yukawa approximation for the special case 
where the stop masses are equal 

In Sect. 6.1.1, we investigated the Yukawa corrections to the triple Higgs vertex and Higgs 
self-energies. For completeness, we give these quantities here in the limit where the stop 
masses in the Yukawa approximation are equal. 

The triple Higgs vertex corrections in this limit are: 



AA 



Yuk 



AA 



Yuk 



A \ Yuk _ A \ Yuk 




AX^uk 



A\Yuk 



0, 



(157) 

(158) 



^Note that the quark field renormahsation factors given in eq. (151) and eq. (152) are sufficient for 
renormahsation purposes but they do not necessarily take into account the full wave function normalisations 
that are required in the most general case to ensure the correct on-shell properties of external quark fields, 
see Ref. [139]. As the quark masses of the final-state fermions are small in the Higgs decays that we are 
considering (i.e., we do not consider Higgs decays into top quarks), in our calculations 6Z^ and SZ^ account 
for the whole wave function corrections. 
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^Ki,. = ,,s:j:3 .5 (159) 



Se^mj I l/i 



2 



AAi^^ = ':!rL .. l-^-244iog^^ (161) 



5 1 m? 



while the Higgs self-energies in this limit are: 

^Yuk - i«^2/l#2 „2 „2 I l/^l ) ' 

f.(i)0i02 ""t / lA-i ""11^ 1 

^Yuk - i«_2/\#2 „2 „2 + o ' ^-^""^^ 




SW^-^^ - 21og(^^^*i^l 1 (164) 



where 



^Yuk — ^Yuk — U, l^iDOJ 



CiiL = Co(0,0,0,m|,m|,M2). (167) 



In this limit, the Yukawa corrections to both the triple Higgs vertex and the Higgs self- 
energies are independent of the CP-violating MSSM phases. 
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